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O \ Abstract 

o 

' We construct exact, entropy satisfying shock wave solutions of the 

■ Einstein equations for a perfect fluid which extend the Oppeheimer- 
^ I Snyder (OS) model to the case of non-zero pressure, inside the Black 

Hole. These solutions put forth a new Cosmological Model in which 
the expanding Friedmann-Robertson- Walker (FRW) universe emerges 
from the Big Bang with a shock wave at the leading edge of the expan- 
^ _ sion, analogous to a classical shock wave explosion. This explosion is 

\^ • large enough to account for the enormous scale on which the galaxies 

, and the background radiation appear uniform. In these models, the 

I shock wave must lie beyond one Hubble length from the FRW center, 

■ this threshhold being the boundary across which the bounded mass 
, lies inside its own Schwarzshild radius, 2M/r > 1, and thus the shock 

wave solution evolves inside a Black Hole. The entropy condition, 
which breaks the time symmetry, implies that the shock wave must 
weaken until it eventually settles down to a zero pressure OS inter- 
I face, bounding a finite total mass, that emerges from the White Hole 

^ ' event horizon of an ambient Schwarzschild spacetime. However, unlike 

shock matching outside a Black Hole, the equation of state p = 
the equation of state at the earliest stage of Big Bang physics, is dis- 
tinguished at the instant of the Big Bang — for this equation of state 



o 



X 



■ alone, the shock wave emerges from the Big Bang at a finite nonzero 



speed, the speed of light, decelerating to a subluminous wave from 
that time onward. These shock wave solutions indicate a new cosmo- 
logical model in which the Big Bang arises from a localized explosion 
occurring inside the Black Hole of an asymptotically flat Schwarzschild 
spacetime. 
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1 Introduction 



In the standard model of cosmology based on a critically expanding, {k = 0), 
Priedmann-Robertson- Walker (FRW) metric, the universe is infinite at each in- 
stant after the Big Bang, 12H l2l]. The Hubble constant, which measures the 
recessional velocity of the galaxies, applies to the entire flat space R'^ — that is, it 
applies to an entire universe of infinite mass and extent — at each fixed positive 
time in the standard model. In this paper we present a new cosmological model 
in which the expansion of the galaxies is a bounded expansion of finite total mass, 
and the Hubble law applies only to a bounded region of spacetime, (c.f. |2B). If 
the observed expansion of the galaxies actually only applies to a localized region 
of spacetime, then it follows that there must be a wave at the leading edge of the 
expansion. Thus to replace the assumption in the standard model that the uni- 
verse is infinite at each instant after the Big Bang, we are led to models in which 
the expansion emerges from an event that is more similar to a classical shock wave 
explosion, than it is to the usual scenario of the Big Bang. 

In previous work, we constructed such models by matching FRW metrics to 
standard Tolman-Oppenheimer-Volkoff (TOV) metrics, (the metric for a static 
fiuid sphere in general relativity), across a shock wave interface, jl4| I21j. In that 
work we derived an upper limit on the distance that a shock wave could be from 
the FRW center, and this distance turned out to be closer than astronomical 
observations suggest it could be. (Astronomical observations currently extend out 
to approximately 1.5 Hubble lengths.) In this paper we begin by showing that 
this upper bound on the shock position that we identified in [2J for standard 
FRW-TOV shock matching, is exactly equal to one Hubble length. Indeed, we 
show that in order for the shock position to lie beyond one Hubble length in an 
FRW metric, it is necessary that the spacetime beyond the shock wave lie inside a 
Black Hole. Thus our previous shock matching limit of one Hubble length results 
from the fact that a standard TOV metric cannot be continued into a Black Hole, 
except in the special case when the pressure is zero, (we proved this in JSI)- With 
this motivation, we here derive a new class of gravitational metrics that we call 
the TOV metric inside the Black Hole. (In contrast to the standard TOV metric, 
the TOV metric inside the Black Hole is dynamical.) Based on this, we develop 
a theory of FRW-TOV shock matching inside the Black Hole, and we use this to 
construct a new class of exact shock wave solutions of the Einstein equations, in 
which a shock wave is incorporated into the FRW metric at distances arbitrarily 
far beyond the Hubble length. This demonstrates that the limit in our previous 
work |21j . that the shock position must lie within one Hubble length of the FRW 
center in FRW-TOV shock matching, can be overcome. 

In the exact solutions constructed in this paper, the expanding FRW universe 
emerges behind a subluminous blast wave that explodes outward from the FRW 
origin f = at the instant of the Big Bang t = 0,^ at a distance beyond one 

■^Hcre, f — R{t)r measures radial arclength distance at each fixed time t in FRW 
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Hubble length'^. The distance of one Hubble length is critical in the FRW metric 
because the total FRW mass M inside radius f satisfies ^ < 1, out to exactly 
one Hubble length. Thus, one Hubble length marks the event horizon of a Black 
Hole in a shock wave model in which the mass M is isolated in an asymptotically 
Schwarzschild metric — the TOV metric in our model. (In contrast to the TOV 
metric outside the Black Hole, the TOV metric inside the Black Hole can be 
continued into an event horizon, c.f. the remarks after Theorem El below.) After 
the Big Bang, the shock wave in our exact solution continues to weaken as it 
expands outward, satisfying the entropy condition for shocks all the way out until 
the Hubble length eventually catches up to the shock wave ^. At this instant 
the shock wave lies at the critical distance of exactly one Hubble length from the 
FRW center. From this time onward, the shock wave can be approximated by a 
zero pressure, k = Oppenheimer-Snyder (OS) interface that emerges from the 
White Hole event horizon of an ambient Schwarzschild metric of finite mass. (The 
entropy condition implies that the TOV density and pressure tend to zero as the 
shock interface approaches the critical distance of one Hubble length.) Thereafter 
the interface continues out to infinity along a geodesic of the Schwarzschild metric 
outside the Black Hole. Thus the OS solution gives the large time asymptotics of 
this new class of shock wave solutions that evolve inside a Black Hole. 

One of the surprises in the analysis is that the equation of state that applies at 
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the earliest stage of Big Bang physics, p = yp, is distinguished by the equations, 
and only for this equation of state does the shock wave emerge from the Big Bang 
at a finite nonzero speed, the speed of light, decelerating to a subluminous wave 
from that time onward.. 

These new shock wave solutions of the Einstein equations confirm the math- 
ematical consistency of an FRW universe of finite extent and non-zero pressure 
expanding outward from behind an entropy satisfying shock wave emerging from 
the origin at subluminous speed beyond one Hubble length at the instant of the 



coordinates, where R{t) denotes the cosmological scale factor, and r is the standard FRW 
radial coordinate, c.f. (|2.4|l . 1)2. 5|l below. 

''The Hubble length c/H{t) = depends on the cosmological scale factor R{t), 

and an easy calculation using the Einstein equations shows that in the FRW spacetime, 
the Hubble length increases with time. Thus more and more galaxies pass inside of the 
threshold distance of one Hubble length and come into view as time evolves. The Hubble 
length c/Hq at present time is estimated to be around 10^° light years. 

^In our cosmological interpretation of the FRW metric, we (loosely) identify the motion 
of the galaxies with the motion of the FRW fluid, which is taken to be a perfect fluid with 
nonzero pressure, co-moving with the FRW metric. We show below that, although the 
shock wave moves outward through the galaxies, (r > 0), and the Hubble length increases 
with time, the number of Hubble lengths from the FRW center to the shock wave, (c.f. 
H5.36|l below), as well as the total mass behind the shock wave, both decrease in time in 
these exact solutions, tending to infinity in backwards time at the instant of the Big Bang. 
This is no contradiction because the FRW pressure p is assumed nonzero, c.f. Corollary^ 
below. 
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Big Bang, a prerequisite for early Big Bang physics. Since the shock wave emerges 
from the Big Bang beyond one Hubble length, it would account for the thermal- 
ization of radiation in a region that is initially well beyond the light cone of an 
observer positioned at the FRW center. Thus our attempt to incorporate a shock 
wave beyond one Hubble length has led to unexpected and interesting connections 
between Big Bang Cosmology and Black Holes. But furthermore, we suggest that 
general relativity pretty much forces the qualitative behavior we see here into any 
reasonable model that relaxes the assumption in the standard model that the ex- 
pansion of the galaxies is of infinite extent at each fixed time. (One could say 
that in these new models, the Copernican Principle is replaced by the principle in 
physics that Nothing Is Infinite.) 

In Section [21 we transform the FRW metric to standard Schwarzschild coordi- 
nates, and use this to discuss the connection between the Hubble length and the 
Schwarzschild radius. In Section |31 we construct the extension of the zero pressure, 
k = 0, OS solution to the interior of a Black Hole by using Eddington-Finkelstein 
coordinates to regularize the event horizon of the Schwarzschild metric, [HI I17j . 
We return to these OS solutions in Section |H1 where we argue that the shock wave 
solutions constructed there continue naturally to the OS solution after the solution 
has relaxed to almost zero pressure. 

In Section 21 we construct the TOV metric inside the Black Hole. We refer 
to this metric as TOV because the metric components depend only on the radial 
coordinate f, but, as in the Schwarzschild metric, the TOV radial coordinate is 
timelike inside the Black Hole. It follows that the mass function M is constant in 
each spacelike slice f = const, of the TOV metric inside the Black Hole. 

In Section [21 we develop the theory of shock matching between a k = FRW 
metric and TOV metrics inside the Black Hole. The shock matching techniques 
introduced in [TBI must be modified inside the Black Hole because the conserva- 
tion constraint used in introduces unphysical (characteristic) solutions. The 
analysis leads to the derivation of a system of differential equations that simulta- 
neously describe the time evolution of the shock position together with the outer 
TOV pressure, and for solutions of these equations, the shock interface must lie 
beyond one Hubble length from the center in the FRW metric. One interesting 
feature of the matching is that the radial coordinate f, a timelike coordinate in 
the TOV metric, is identified through shock matching with the FRW spacelike 
coordinate f that measures radial arclength distance from the center at each fixed 
value of the (standard) time t in the FRW metric. Another interesting feature is 
that the mass function M, which is continuous across shocks as a consequence of 
shock matching, has the physical interpretation as a total mass inside radius f in 
the FRW metric, but M has no such interpretation in the TOV metric by itself, 
and in fact we know of no general physical interpretation of the mass function M 
inside the Black Hole. 

In Section[Hlwe formulate the entropy condition, and construct a class of global, 
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entropy satisfying shock wave solutions of the Einstein equations under the simph- 
fying assumption that the FRW sound speed ^ is constant, (that is, we assume 
the FRW equation of state -p = cr/o, a = const. > 0, c.f. JU). This includes 
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the important case a = usually assumed at the earliest stage of Big Bang 
physics. (This is the equation of state in the extreme relativistic limit of free 
particles, and for pure radiation, ^S]-) Remarkably, under a change of variables, 
the shock matching equations of Section [3 reduce to a planar autonomous system 
when a = const., and this system is amenable to global analysis, a requirement 
for the construction of solutions in the large. It is very interesting that the special 
value a = c^/3 plays distinguished role, and at this unique value, the shock wave 
is everywhere subluminous after the Big Bang, but emerges from the Big Bang 
at exactly the speed of light. For all other values of a, we prove that the shock 
speed at the instant of the Big Bang is either zero, or infinite, and is everywhere 
subluminous if and only if o" < c^/3. This is surprising because the equation of 
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state p = yP played no special role in shock matching outside the Black Hole. 

The class of exact solutions in Section describes the global dynamics of an 
FRW universe of finite extent which explodes outward behind an entropy satisfying, 
subluminous shock wave, emerging from the origin, beyond one Hubble length, at 
the instant of the Big Bang. Because the TOV metric inside the Black Hole has 
nonzero density and pressure, it follows that, unlike the OS solutions of Section |21 
these new exact shock wave solutions do not require any part of the empty space 
Schwarzschild solution inside the Black Hole for their construction. On the other 
hand, the fact that the OS solution gives the qualitative large time behavior of 
the solutions independent of a, (a consequence of the entropy condition alone), 
suggests that the qualitative features of these solutions may be generic for a large 
class of equations of state. 

In Section [7| we obtain estimates for the shock position. The conclusion of 
Theorem ini is this: Let to be the first time at which the shock becomes visible at 
the FRW center. In our exact solutions, we assume the FRW equation of state 
p = ap, a = const, < a < 1/3. (We use the convention that we take c = 1 when 
convenient.) For these solutions, the shock wave will first become visible at the 
center f = of the FRW spacetime at FRW time t = to, at the moment when the 
Hubble length satisfies 



H{to) 



7, (1-1) 



where r^, is the FRW position of the shock at the instant of the Big Bang. (The 
actual arclength distance from the center to the shock at the Big Bang will be 
f = 0. ) At this time, the number of Hubble lengths ^/N q from the FRW center 
to the shock wave at time t = to satisfies the bounds 

-1 + 30-- l + 3o- 
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Thus, in particular, the shock wave will still lie beyond one Hubble length at the 
FRW time to when it first becomes visible. Furthermore, the time tcrit > at 
which the shock wave will emerge from the White Hole event horizon, given that 
Iq is the first instant at which the shock becomes visible at the FRW center, can 
be estimated by 



e'i < < e 1+-^ , 

1 + 3o- ~ ^0 ~ 1 + 3o- 

and by the better estimate 

^ tcrit , 3 
e 4 < < 62 , 

" to " 

in the case a = 1/3. For example, (fTHfl . imply that at the OS limit a = 0, 




and in the limit a = 1/3, 

l<V^<l-95, 1.8< — <4.5. 

to 

We conclude in these shock wave cosmological models, that at the moment t=K 
when the shock wave first becomes visible at the FRW center, it must lie within 
4.5 Hubble lengths of the center. Throughout the expansion up until this time, 
the expanding universe must lie entirely within a Black Hole — the universe will 
eventually emerge from this Black Hole, but not until some later time tcrit: where 
tcrit does not exceed 3^=,,. 

Concluding remarks are made in the final section. 



2 The Hubble Distance and the Schwarzschild 
Radius 

According to Einstein's Theory of General Relativity, all properties of the gravita- 
tional field are determined by a Lorentzian spacetime metric tensor g, whose line 
element in a given coordinate system x = (x", ...,x^) is given by 

ds"^ = Qijdx'dxK (2.1) 

(We use the Einstein summation convention whereby all repeated up-down indices 
are assumed summed from to 3.) The components gij of the gravitational metric 
g satisfy the Einstein equations, 

G^^ = kT'^, = {pc^ + p)w'w^ + pg'\ (2.2) 
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where we assume the stress-energy tensor T of a perfect fluid. Here G is the 
Einstein curvature tensor, 



K = (2.3) 

is the couphng constant, Q is Newton's gravitational constant, c is the speed of 
light, pc^ is the energy density, p is the pressure, and w = {w^ , ...,w^) are the 
components of the 4- velocity of the fluid, c.f. |21]. We use the convention that we 
take c = 1 and Q = 1 when convenient. 

In this section we consider the FRW metric, a spacetime metric whose line 
element takes the form 

^ _^^2 ^ ^(^)2 (^-^ + r'dn"^^ , (2.4) 

where R{t) is the so-called cosmological scale factor, and dQ'^ = dO^ + siii?9d(fP' 
denotes the line element on the unit 2-sphere. The FRW metric describes the time 
evolution of a three dimensional space of constant scalar curvature, (the i=const. 
surfaces), and k can be assumed to take one of the values {—1, 0, 1} via a rescaling 
of the radial coordinate r. Radial distance at each fixed time in the FRW metric 
is measured by f, where 

f = R{t)r, (2.5) 

and it is standard to rescale time so that R goes from to 1, where R = 0, t = 
corresponds to the Big Bang, and R = 1, t = to corresponds to present time. 

The Hubble constant H, (which actually depends on t), is given in terms of 
-R(t) by 

H . H(t) = f||. (2.6) 

and we let Hq = i?(to) denote the present value of the Hubble constant. 

The Hubble length c/H, (the reciprocal of the Hubble constant when we take 
c = 1), gives a length scale determined by the expansion rate of the galaxies, and 
can be interpreted as the travel distance for a light ray starting at the Big Bang 
and ending at time t. (That is, the age of the universe is approximately H^^ ~ 10^*^ 
years, so the Hubble length c/Hq « 10^° light years, is a measure of the distance 
light travels during this time interval, The Hubble length is thus a measure 

of the distance to the furthest objects that can be seen in the universe at a given 
time. (Estimates for this distance in terms of radial arclength f at fixed FRW time 
t are given in Theorems ^ and |21 below.) 

Putting the metric ansatz (|2.4j) into the Einstein equations (|2.2)) gives the 
equations for the FRW metric, ^24,, 
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) 



2 



K k 



3^ ii2 



(2.7) 



and 



P = - 



3{p + p)H. 



(2.8) 



The unknowns it!, p and p are assumed to be functions of the FRW coordinate 
time t alone, and "dot" denotes differentiation with respect to t. In this paper, we 
focus on the case of critical expansion, the case /c = 0, [Tl [^. 

When k = 0, the t = const, surfaces in the FRW metric are infinite, flat, 
Euclidean 3-space. It follows that the standard model of cosmology based on a 
critically expanding FRW metric, (Ij, implicitly assumes that the expansion rate 
H of the universe is constant all the way out to infinity at each fixed time t in the 
FRW metric — this even though, based on lookback time, we can in principle only 
observe the universe out to about 10^° hght years. 

If, on the other hand, there is a shock wave at the leading edge of the expansion 
of the galaxies, then the following question presents itself: What is the critical 
radius fcrit at each fixed time t in a = FRW metric such that the total mass 
inside a shock wave positioned beyond that radius puts the universe inside a Black 
Hole? Indeed, there must always be such a critical radius because the total mass 
M{f,t) inside radius f in the FRW metric at fixed time t increases like f^, and 
so at each fixed time t we must have f > 2M(f, t) for small enough f, while 
the reverse inequality holds for large f. Thus for every time t there must exist a 
smallest f = fcrit for which fcrit = '^M{fcrit,'t). This critical radius then marks the 
Schwarzschild radius in the metric that lies beyond the shock wave when it is at 
that position. We will presently show that when k = 0, fcrit equals the Hubble 
length, and this explains the shock matching limit uncovered in our previous work 
in terms of the Hubble length, (c.f. Section 6, 125 )■ That is, we cannot match 
a critically expanding FRW metric to a classical TOV metric beyond one Hubble 
length without continuing the TOV solution into a Black Hole — and we showed in 
|15j that the standard TOV metric cannot be continued into a Black Hole. 

At this point we make a note on terminology: we say that a radial solution of 
the Einstein equations is inside the Black Hole provided that ^M- > 1. If further, 
the radial coordinate f is always increasing along timelike curves, then we also 
refer to such a region ^M- > 1 as a White Hole, but we will nevertheless use the 
term Black Hole whenever ^M- > 1.. For example, in the Kruskal development of a 
Schwarzschild Black Hole there are four quadrants, and one quadrant contains the 
Black Hole singularity and the opposite quadrant, (its time reversal), contains the 
White Hole singularity, 9 . Both quadrants lie within the region where ^y- > 1, 
and the quadrant with the Black [resp. White] Hole singularity at f = has the 
property that all timelike geodesies end [resp. begin] at the singularity in finite 
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proper time^. Thus if fcrit = c/H, and the mass function M is continuous across 
shock waves, (both of which we demonstrate below), then it foUows that one can 
incorporate a shock wave into an FRW metric beyond one Hubble length only if 
the metric beyond the shock wave lies inside a Black Hole, where ^^^i'''*^ > 1. 

We now verify that fcrit = c/H, and that the Hubble length is the limit for 
FRW-TOV shock matching outside a Black Hole, c.f. |2J. To this end, write 
the FRW metric (|2.4j) in standard Schwarzschild coordinates R = {f,t) where the 
metric takes the form 

ds^ = -B{r, t)df + A{r, t)-^df'^ + r'^d^'^, (2.9) 
and the mass function M(f, t) is defined through the relation 

A = l-^. (2.10) 

r 

(We have set Q = c = 1, and our notation is to denote the standard, non-angular 
FRW coordinates by x = {x^,x^) = {t,r), but TOV coordinates are denoted by 
because we work inside the Black Hole where ^ < and r 
is the timelike coordinate.) It is well known that a general spherically symmetric 
metric can be transformed to the form 1)2. 9|) by coordinate transformation, |24j . 
To obtain A and B for the FRW metric set f = Rr, so that 

, df f ■ 
Using this in the FRW line element we find 

ds'^ = -Cdt^ + Ddf^ + 2Edfdt + f^dl^^ (2.11) 

where 



C = '-"f-f", (2,12) 



1 

1 ~kr'^ 
Hf 



D = (2.13) 



^By itself, Einstein's theory is time reversible, the direction of forward time is left unde- 
termined, and extra conditions like entropy conditions for shocks are needed to determine 
the direction of forward time. In the region 2M/f > 1 outside the event horizon, the 
Schwarzschild metric is symmetric under time reversal in the sense that the time reversal 
of a geodesic remains a geodesic. However in the region 2M/f < 1, the time reversal of 
the Schwarzschild metric takes the Black Hole region beyond t = +oo to the White Hole 
region before t = —oo. Thus without extra conditions, the expanding White Hole is just 
as admissible a solution of the Einstein equations as is its time reversal, the collapsing 
Black Hole. 
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Now define the time coordinate t = t{t, r) by 

dt = {ipC)dt - {^pE)dr, (2.15) 
where V' is an integrating factor to be determined by the equation 

^(V^C7) + |(V.i?)=0. (2.16) 

Equation 1)2. 16(1 imphes that dt is an exact differential, and the choice of E and C 
applied to (|2.15|) removes the cross term from the metric (|2.11|) . In (f, t) coordi- 
nates the FRW metric takes the form 

= --^dP +(^D + ^^ df^ + r^d^\ (2.17) 
Using (l2T2|l - (l2ni1) in (tOTl) gives 

1 I l-kr"^ 1 ^ r 1 I .^2 , ^2 ,o2 



ds' = dt' + | i_,,2_^2,2 ^ + ^'dn^- (2-18) 

Comparing 1)2. 18() with ()2.9() we obtain 




B = — I . . . ,_„.2 ■ (2-20) 



1 — kr"^ — (rH) 
Note that the characteristic curves for ()2.16p are given by 

dr _C _ -l + fcr2 + i/2^2 
dt ~ E~ ' 



(2.21) 



the RHS being a smooth function for all t > for which H ^ ^ f. (By (|2.7j) . 
the condition H = ^ ^ holds for all t > when k = 0,-1, and it holds except 
at the turning point = in the case when k = +1, c.f. 24 .) It follows that a 
solution ip of ()2.16|) exists in a neighborhood of any point. In fact, we have more: 
we can assign arbitrary initial values for -0 on any surface that is non-characteristic 
for (|2.1H|) . and then use these values to solve (|2.16|) for -0 in a neighborhood of 
that surface, c.f. jl4l I24j . By (|2.2U|) . this implies that we can assign arbitrary 
(non-negative) values for B on any non-characteristic surface. We conclude that 
(|2.18() - (|2.2U() gives the FRW metric in standard Schwarzschild form, the coefficient 
of df2 is determined independently of and values of B can be freely assigned 
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(locally) on any smooth non-characteristic surface in (t, r) space, thus determining 
initial values for equation (|2.16j) . 

Moreover, using ((TT^ in (HHUl), it follows that 



= I r p{t)s^ds. (2.22) 
2 Jo 

Since in the FRW metric f = Rr measures arclength along radial geodesies at fixed 
time, we see from (|2.22p that M(i, f) has the physical interpretation as the total 
mass inside radius f at time t in the FRW metric. 

From here on, we restrict to the case of critical expansion, the case fc = 0. In 
this case, since ^ = 1 — A, one sees from (|2.19j) that f = is equivalent to 
'^M- = 1, and so the following equivalences are valid at any fixed time t: 



f = H-^ iff ^ = 1 iff ^ = 0. (2.23) 

r 

We conclude that f = is the critical length scale for the FRW metric at fixed 
time t in the sense that "^M- — 1 changes sign at f = , and so the universe lies 
inside a Black Hole beyond r = H~^, as claimed above. 

From the FRW equation 1)2. 7|) . it also follows that f = H^^ iff 

This explains the condition 

9 3 1 , , 

given in equation (6.2) of !21^ for the maximal position of a shock wave in FRW- 
TOV matching. This condition was derived in [2XJ from a different point of view. 
Thus (|2.23j) . (|2.24)) provides a physical interpretation for the bound (6.2) in j21j : 
namely, that the shock wave lies inside the Hubble length. 

The Hubble length rcrit — Ji '^^ ^Iso the critical distance at which the outward 
expansion of the FRW metric exactly cancels the inward advance of a radial light 
ray impinging on an observer positioned at the origin of a /c = FRW metric. 
Indeed, by ()2.4|) . a light ray traveling radially inward toward the center of an 
FRW coordinate system satisfies, 

c^dt^ = R^dr'^, (2.26) 

so that 
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— = Rr + Rr = Hf - c = H{f - — ) > 



(2.27) 



if and only if 



c 



Thus the arclength distance from the origin to an inward moving Ught ray at fixed 
time t in a A; = FRW metric will actually increase as long as the light ray lies 
beyond the Hubble length. An inward moving light ray will, however, eventually 
cross the Hubble length and reach the origin in finite proper time, due to the 
increase in the Hubble length with time. 

We now calculate the infinite redshift limit in terms of the Hubble length. It 
is well known that light emitted at (tejJ'e) at wavelength Ag in an FRW spacetime 
will be observed at {tQ^ro) at wavelength Aq if 



Moreover, the redshift factor z is defined by 



Thus, infinite redshifting occurs in the limit i?e — ^ 0, where i? = 0, t = is the 
Big Bang. 

Consider now a light ray emitted at the instant of the Big Bang, and observed 
at the FRW origin at present time t = Iq. Let r^o denote the FRW coordinate at 
time i — > of the furthestmost objects that can be observed at the FRW origin 
before time t = to. Then r^o marks the position of objects at time t = whose 
radiation would be observed as infinitly redshifted, (assuming no scattering). Note 
then that a shock wave emanating from f = at the instant of the Big Bang, will 
be observed at the FRW origin before present time t = to only if its position r at 
the instant of the Big Bang satisfies r < Too- We now estimate . 

First, from (|2.26|1 it follows that an incoming radial light ray in an FRW metric 
follows a lightlike trajectory r = r{t) if 





and thus 




(2.28) 



We now prove the following theorem: 
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Theorem 1 If the pressure p satisfies the bounds 

0<p<^p, (2.29) 

then for any equation of state, the age of the universe to cind the infinite red shift 
limit Too are bounded in terms of the Hubble length by 

2^ ^ *° S ^l^' 

12 , , 

TT < < TT- (2-31) 

(We have assumed that R = when t = and -R = 1 when t = to and H = Hq. 
Note that p = ^/o is the extreme relativistic hmit of free particles, as weh as the 
equation of state for pure radiation, ^^.) 

Proof: Integrating (|2.7|1 we obtain 

R{t) = e~-lt v-"'*, (2.32) 

and, from (|2.8|) we have 

dp 



Thus, if we know p as a function of p, then we can integrate ()2.33() to obtain p as 



a function of t, and then use this in (|2.32j) . which can then be used in (|2.28|) to 
calculate Too- Assuming (|2.29|) . we can estimate (|2.33p by 



which leads to 



AKt"^ - ^ - 3Kt^ ■ 
Since (ETj) is H"^ = §p when A; = 0, (lO^ gives 



3 4 , , 

<P<— ^. (2.35) 



1 2 

— < H < —, (2.36) 
2t - - 3t' ^ ^ 



which implies H2.3U() . To estimate R{t), write ()2.36|) in the form 



^ < ^InR < — , (2.37) 
2t - dt - 3t ^ ' 



which integrates to 
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Finally, using ((T^ in ^TM gives 

(t) *^'-^/o (t) ^'■''> 

which leads directly to ()2.3ip . □ 

The next theorem gives closed form solutions of the FRW equations (|2.7j) , (|2.8j) 
in the case when the sound speed ^/a^ = constant. These solutions are the starting 
point of the exact shock wave solutions constructed in Section 6. As a special case 
we recover the bounds in (|2.3()|) and (|2.31|) from the cases a = and 1/3. 

Theorem 2 Assume k = and the equation of state 

p = ap, (2.40) 

where a is assumed constant, 

< a < 1. 

Then, (assuming an expanding universe R> 0), the solution of system \2. 7| ), \2. t^) 
satisfying R = at t = and R = 1 at t = to is given by, 

P=M^W^ (2-41) 

W, = 'f- (2-43) 

Moreover, the age of the universe io (^n-d the infinite red shift limit r^o are given 
exactly in terms of the Hubble length by 



From 1)2. 45() we conclude that a shock wave will be observed at the FRW origin 
before present time t = to only if its position r at the instant of the Big Bang 
satisfies r < jq^-^- Note that Too ranges from one half to two Hubble lengths 
as cr ranges from 1 to 0, taking the intermediate value of one Hubble length at 
a = 1/3, c.f. (lOTTl . 

Proof: Formulas (|TlT|) - (|nH|) follow directly from and agree with 

the formulas given in p2|- Differentiating ()2.42|) at t = to gives ()2.44() . and using 
l^rm and (ITTO in fPTM gives flJ^ . □ 
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Corollary 1 If p = ap, a = const. > 0, then the total mass inside radius 
r = const, (that is, inside a ball whose boundary is comoving with the galaxies), 
decreases in time. 

Proof: Using ([2lTI) - ([02|l in ((2221), it follows that 



M 



- p{t)s'ds 



2r^ 



9(1 + 



-2a- 



(2.46) 



so < if cj > 0, 



3 The OS Solution Inside the Black Hole 

The simplest model of a localized FRW metric contained within a shock boundary 
that lies beyond the Hubble distance, is one in which p = 0, and the FRW metric is 
matched to the Schwarzschild metric at a contact discontinuity interface positioned 
inside the Black Hole of the Schwarzschild metric. This poses the problem of 
extending the OS solution smoothly into the interior of a Schwarzschild Black 
Hole. In this section we construct the OS solution inside the Black Hole in the 
case /c = 0, making the point of discussing it from the FRW point of view, in terms 
of the Hubble length. (See ITOl and l9l for a discussion of the case A; < 0.) The 
k = Q OS solution inside the Black Hole will give us a point of comparison for the 
p 7^ shock wave models constructed in Sections 5 and 6, and we argue in Section 
7 that the (expanding) OS solution also describes the large time asymptotics of 
these shock wave models. When p 7^ 0, the resulting interface is a shock wave, 
and an entropy condition breaks the time symmetry. In our examples of Section 
6, it is the outward expanding solution (explosion) that globally meets the entropy 
condition, (not the inward collapse), thus making the cosmological interpretation of 
the model relevant. In contrast, the interface in the OS solution is time reversible, 
but the collapsing solution is the one that is relevant in the standard interpretation 
of the OS solution as a pressureless sphere collapsing into a Black Hole, [TUlin]. 

Our method is to match the fc = FRW metric to the empty space Schwarzschild 
metric written in Eddington-Finkelstein (EF) coordinates P|, across a shock in- 
terface'' that lies beyond the Hubble length on the FRW side of the shock. As in 
a classical explosion, we assume that the FRW metric lies inside a bounded region 
behind an outgoing shock interface, and for the OS solution we assume that the 
Schwarzschild metric describes the spacetime beyond the interface. Thus in the 
OS model, the shock wave (contact discontinuity) marks the leading edge of the 
expansion of the FRW metric. 

^The interface in the OS solutions is, in the language of gas dynamics, a contact 
discontinuity, which is time reversible because neither mass nor momentum cross the 
interface. 
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The EF coordinates provide a regularization of the Schwarzschild metric at 
the event horizon, and the mapping from the EF spacetime to the Schwarzschild 
spacetime is a 1-1 mapping, [5]. Both coordinate systems cover the region of the 
Schwarzschild metric outside the event horizon, together with that portion of the 
Black Hole that is coordinatized by the Schwarzschild coordinates. Here, since we 
are dealing with an expanding FRW solution, we work with the time reversal of 
these metrics, which covers the region beyond the event horizon, together with 
the White Hole region inside the Black Hole We show that with this latter time 
orientation, we can match the EF metric to the spacetime described by a critically 
expanding FRW metric across a discontinuity at finite radius, and its time reversal 
provides the corresponding matching to the critically contracting FRW metric. 

We return to these OS solutions in Section 6 where we argue that the shock 
wave solutions constructed there continue naturally to the OS solution after the 
TOV density and pressure tend to zero, assuming that the shock has relaxed to 
a sufficiently weak wave. But we find this OS model interesting in its own right 
because it not only provides the simplest model in which the expansion of the 
galaxies corresponds to the expansion of a finite total mass with a wave at the 
leading edge of the expansion, but it also embeds the Big Bang singularity of an 
FRW metric within the singularity of a larger spacetime — the larger spacetime 
being the empty space Schwarzschild solution, and its singularity being the White 
Hole singularity that lies inside the event horizon in the Kruskal development of the 
Schwarzschild metric. Such an embedding is possible only under the assumption 
that the Hubble law applies to a bounded region of spacetime, because the infinite 
FRW metric cannot be matched to the Schwarzschild metric. 

In EF coordinates, the Schwarzschild line element takes the form 



ds^ = -AdP - 2dfdi + f^dl^^ (3.1) 

where 



2M , , 

A = l-—, 3.2 
r 

and M is the constant mass of the Schwarzschild solution. Indeed, using (|2.1H) - 
()2.16|) . the time coordinate i that eliminates the mixed term in 1)3. 1|) satisfies 



dt = ^jAdi - ipf, (3.3) 

where ■i/' satisfies 



|V' + ^(V^) = 0. (3.4) 
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The solution to (jS31) is 



iP = IM, (3.5) 



t = t -f -ln\f -2GM\, (3.6) 
and in (f, f)-coordinates the metric takes the standard Schwarzschild form 

ds^ = -AdP + \dr^ + r^dn\ (3.7) 
A 

Thus 1)3. 6|) defines a 1 — 1 mapping {t, r) — > (f, t) taking {f > 2M} x R to itself, 
and another 1 — 1 mapping (t, f) — > (f, t) taking {0 < f < 2M} x R to itself. 
This verifies the claim that the EF metric covers two of the four quadrants of the 
(maximal) Kruskal development of the Schwarzschild metric in a 1 — 1 fashion. 
Unlike the Schwarzschild metric, the EF metric is smooth at f = 2M, and this is 
reflected in the singularity of the transformation (|3.6j) at f = 2M. Note that in 
(f, t) coordinates, the vector field ^ is timelike, spacelike, when f < 2M, f > 2M, 
respectively; while in (i, f) coordinates, the vector field ^ is lightlike. This is 
no contradiction because ^ points along the level curves of its complementary 
coordinates in a given coordinate system. On the other hand, the vector fields ^ 
and ^ both point along the level curves of f, and hence retain the same character 
as spacelike or timelike according to the sign of A. 
The A; = FRW metric is given in (|2.4j) as 



ds2 ^ _^^2 ^ ji^tf^r"^ + R{tfr^dn^, (3.8) 

where R{t) is the cosmological scale factor. To obtain a matching of the two 
metrics (|3.1|) . (|3.8|) . that is smooth at the Schwarzschild radius, we match them in 
{t, f) coordinates. To start, first match the spheres of symmetry by setting, (c.f. 

mm), 

f = Rr. 

Writing (|2.4|) in (t, f) coordinates yields the following form of the FRW metric 
derived in 1)2. llj) - 1)2. 14(1 above upon setting k = 0: 



ds"^ = -dt^ + \dr^ - 2Hfdfdt + ^r'^dt^ | + r^d^' 



(3.9) 
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We now construct the transformation t = i{t, r) by finding functions b{t, r) and 
(j){t,f) such that b and cp satisfy 



di = (j)dt + (j)bdr, 



(3.10) 



and 



0. 



(3.11) 



Equation ()3.11|) imphes that dt is an exact differential, and thus (|3.1U|) defines the 
coordinate transformation for t as a function of t and f. In order to find (j) and b, 
we start from (|3.1()|) and write 



dt = i |dt - 06df| , 



and so 



(3.12) 



Using these in (|3.9j) we obtain 



(3.13) 



ds^ = J_(4f2-iUt' + -<!-M-i + ^r 



62 I i?2 



1 + 6 





Hf } dtdr + 



(3.14) 



The outline of our procedure for matching the components of the metric (|3.14j) to 
the components of the metric 1)3. 1() is as fohows: We first set the coefficient of df2 
.141) and solve for the function b. We then match the coefficients 



equal to zero in 

of the cross terms dtdf by setting the dtdf coefficient in ()3.14|) equal to —2 and 
then solve for the value of cp on the shock surface. Using the values of b and (p 
in the matching of the dP coefficients, we obtain a formula for the shock surface. 
Finally, we show that the shock surface is non-characteristic for the PDE (|3.11() . 
implying that we can use the initial condition for (p to solve for (p in a neighborhood 
of the shock surface, and thereby determine a nonsingular coordinate t = t{t,r). 
To begin, set the coefficient of df"^ equal to zero in (|3.14|) to obtain 



1 + 6 {6 (-1 + H^f^^ + 2Hr} = 



0. 



(3.15) 
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This simplifies to 



{Hrb+iy = 6^ 

so that 



^ (3.16) 



Hf+l 

In order to ensure that b, given in H3.16() . is nonsingular and of a single sign, we 
choose the plus sign if i? > 0, (the case of a White Hole explosion), and the minus 
sign if i? < 0, (the case of collapse to a Black Hole). Matching the cross terms in 
JSH) and yeilds 

-{-6(-l + -H'2f2) -/7f} = -2. (3.17) 



Using (|3.15j) . we obtain 



Hr + ^=^1 (3.18) 



at the shock surface. Here, consistent with (|3.16)) . we take +1 if H = ^ > 0, and 

R 

R 



-1 if H = ^ < 0. Matching the dP terms in (|XT|) and (|XTH) gives 



-A = ^{H'r'-l). (3.19) 
Since A; = and p = 0, the FRW equations (fT7|) . (|^ imply 



= ^p. (3.20) 

Using (|8.20j) and (|3.2|) together with (/)^ = 1 in (|3.19|) yields the equation for the 
shock surface, 



M = ^p{t)r\ (3.21) 

Solving (|3.2H) for f as a function of t gives the shock surface f = f{t) in the FRW 
coordinate. It follows from the FRW equations 1)2. 7|) . 1)2. 8|) in the case p = 0, that 
R{t)'^p{t) is constant on FRW solutions, [2^. Using this in (|3.21|1 implies that 
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r = ro = const, describes the interface in the OS solution. For example, if we 
choose R{to) = 1 at present time t = to, then 



r 6M 1 

is constant on the interface, and thus f(t) = R{t)rQ describes the motion of the 
interface in the OS solution. The fact that r = const, along the interface implies 
that the interface moves with the galaxies of the p = FRW metric, and hence nei- 
ther mass nor momentum crosses the interface. From this we can also see that the 
shock surface is non-characteristic for the PDE (|3.11j) . Indeed, the characteristic 
curves for (|3.1H) satisfy 

on the shock surface, where the 1 agrees with the sign of R. Thus we see that 
^ = {Rro)tl along the characteristic, and ^ = Rro along the shock surface, so the 
two can never be equal. We conclude that the shock surface is never characteristic 
for equation (|3.1H) . and hence (j), as well as the coordinate t, can be defined in a 
neighborhood of the interface. 

Note that the shock surface equation (|3.2H) tells us that f goes from zero to 
infinity as the FRW density goes from infinity to zero. For example in the case 
R > 0, this tells us that the FRW universe starts at the White Hole singularity 
of the EF metric at the instant t = of the Big Bang. From this time onward, 

the shock surface f(t) = (^^^^) continues out until the the Hubble distance 

1 /2 

H^^ = (^^j^l^y) ' ('^■^^ (Eiini))) catches up to the shock surface on the FRW side 
of the shock, which happens at the FRW time tg when f{ts) = 2M, i.e., when the 
shock surface lies on the event horizon of the outer EF metric. As the FRW time 
increases from tg to oo, the interface continues on out to infinity, staying inside 
the Hubble length on the FRW side of the shock, and outside the Black Hole on 
the EF side. 

Note also that the coordinate f is a spacelike coordinate that measures arc- 
length distance in the FRW metric at fixed time t in the FRW coordinate system 
(t, f) on the FRW side of the shock, but on the EF side of the shock, the coordinate 
f is lightlike in EF coordinates {t,f), and changes from timelike to spacelike in 
standard Schwarzschild coordinates (f , t) as the shock surface passes through the 
event horizon of the EF metric — even so, the coordinate identification (t, r) — > 
(t, f) is regular in a neighborhood of the shock surface for all shock positions 
< f < cxD. This is no contradiction, because the coordinate vector field ^ 
remains undetermined until a choice of complementary coordinate is specified. 

In conclusion, we have that the interface between the FRW and Schwarzschild 
metrics in the A; = OS solution is a contact discontinuity that traverses a geodesic 
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of the Schwarzschild metric. The total mass of the FRW metric behind the shock 
wave interface is finite and constant in time, and emerges from the White Hole 
singularity of the ambient Schwarzschild metric at f = 0, the instant of the Big 
Bang. After a finite proper time, the solution continues out through the event 
horizon at Schwarzschild time t = —oo, and continues to expand forever into the 
asymptotically flat Schwarzshild spacetime outside the Black Hole, where ^M- < 1. 
In this regime the solution agrees with the A; = OS solution given in jl7p. 

4 The TOV Solution Inside the Black Hole 

The standard Tolman-Oppenheimer-Volkoff (TOV) metric models a static fluid 
sphere in general relativity. It was proved in that the standard TOV metric 
for a static fluid sphere, by itself, does not admit Black Holes, and can only exist 
when ^M- < 1. In this section we derive the equations for the analogue of the TOV 
solution inside the Black Hole. 

The TOV solutions inside the Black Hole are used in Sections|Sland|Hlto extend 
the OS solution of Section IHl to the case of non-zero pressure, inside the Black 
Hole. When p ^ 0, energy and momentum cross the interface, and so the contact 
discontinuity of the OS solution must be replaced by a shock wave discontinuity. 
Thus in order to extend our shock matching techniques beyond the Hubble length 
when p ^ 0, we must replace the outer Schwarzschild metric of the OS solution 
with a metric that contains matter, and satisfies 2M/f > 1. The TOV metric 
inside the Black Hole is the simplest metric that satisfies these conditions. We 
refer to it as TOV because the components depend only on the radial coordinate 
f, as in the standard TOV metric, but now f is timelike. (As in the Schwarzschild 
metric, the roles of space and time are interchanged inside the Black Hole.) In 
shock matching with FRW metrics, this new TOV metric inside the Black Hole 
will play the role of a transitional solution that mediates the mass flux across the 
shock interface during the time after the Big Bang when the densities are large, 

®It is well known that the spacetime coordinatized by either the EF or Schwarzschild 
coordinates is not a geodesically complete spacetime. That is, from the point of view of 
the (complete) Kruskal development of the Schwarzschild spacetime, the EF metric covers 
only two of the four quadrants determined by the event horizon in the Kruskal diagram; 
namely, the standard quadrant outside the Black Hole together with the quadrant which 
contains the Black Hole singularity-or else under time reversal, the standard quadrant 
outside the Black Hole together with the region containing the White Hole singularity. 
On the other hand, the OS solution is a geodesically complete spacetime because, from 
the Kruskal point of view, it consists of the Schwarzschild solution on one side, and the 
FRW solution on the other side, of the geodesic defined by the trajectory of the interface 
of discontinuity. Thus, by removing a neighborhood of either the White Hole or the Black 
Hole event horizon, the OS matching also eliminates the incomplete geodesies that emerge 
from t = — oo, respectively t = +oo,, depending on the time orientation of the solution, 

i- 
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and up until the time when the solution has settled down to a zero pressure OS 
expansion. 

The usual ansatz for a TOV metric takes the form 

ds^ = -Bdf + A-^dr^ + f^dn^ (4.1) 
where A(f) and -B(f) depend only on the coordinate f. Here 

^(f) = l-^. (4.2) 

r 

and when f > 2M, M{f) has the interpretation as the total mass inside the ball 
of radius f. Thus our assumption that the TOV metric lies inside the Black Hole 
or inside the Schwarzschild radius is equivalent to the assumption that 

A<0. 

We now obtain the Einstein equations for a perfect fluid under the assumption 
that the fluid is co-moving with respect to a metric of form (|4.1|) . assuming that 
A and B depend only on f, but now assuming that ^ < 0. The stress tensor for a 
perfect fluid takes the form, 



Tij = {p + p)wiWj + pgij, (4.3) 

and since f is the timelike coordinate for the TOV metric in (f , f)-coordinates when 
A < 0, the assumption that the fluid is co- moving with the TOV metric inside the 
Black Hole implies that, c.f. |24j . 




In this case we obtain 



00 



T— 



(P+P) 



-A 



+ P 



P_ 

A' 



(4.5) 



and 



Tn = Trt= pgn = -pB. (4.6) 

Using MAPLE we find, 

TAB' -B + AB 

(-^00 = <J^ff — r-TB ' l^.^j 

rAB 
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Gi^ = Grt = -^{rA' -l + A). (4.8) 
From (|4.5|1 - (|4.8|) . the Einstein equations Gqo = hTqq and Gu = kTu reduce to 

B' 1 /1-A\ p ,^ ^ 

and 

\ — A 

A' = — z hKpf, (4.10) 

r 

respectively.^ From equation (5.45), we also have 

-P' = (p + P)^ln{^}, (4.11) 

which simphfies to 

Using (|4.2|) . we obtain the system 



p + pA' p + p(l-A __\ 

P = ^:4=^l^ + "^7' ^'-''^ 

1 — A 

A' = — z \-Kpr, (4.14) 

r 

B' I /1-A\ p , , 

Alternatively, using the unknown = 1 — ^ instead of A, (a variable convenient 
for our subsequent analysis), we obtain the equivalent system, 

= i±PJ^, (4.16) 

N' = -\^ + Kpf\, (4.17) 



J. 

B' I (N 



B N -I 



Kp\. (4.18) 



^Beware that MAPLE's convention is that ~G maple — G, so the Einstein equations 
are -Gmaple = kT. 
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Note that the essential reason why the equations for the TOV metric inside the 
Black Hole take a different form than the standard TOV equations "outside the 
Black Hole" , is that the assumption of co-moving puts the nonzero component of 
w on the timelike coordinate f when ^ < 0, and on the timelike coordinate i when 
^ > 0. We conclude that system (|4.1H|) - (|4.18|) defines a new class of gravitational 
metrics which describe spacetimes that evolve inside the Black Hole. 

The TOV equations inside the Black Hole describe a time dependent metric 
in which the metric components, together with the fluid variables, are constant at 
each fixed value of the TOV timelike coordinate f. Like the FRW metric, the TOV 
metric inside the Black Hole describes a fluid with pressure that emerges from a 
White Hole singularity at f = 0, the instant of the Big Bang, with one important 
difference: unlike the FRW metric, for the TOV metric inside the Black Hole, the 
total mass is constant on each spacelike slice f = constant. In the next section we 
match the TOV metric inside the Black Hole to k = FRW metrics across shock 
interfaces in order to make the expanding FRW universe finite, in the sense that 
the total mass on each spacelike hypersurface is finite. 

5 Shock Matching Inside the Black Hole 

In this section we derive the equations that describe the matching of a general 
k = FRW metric to a TOV metric inside the Black Hole, at a shock wave 
interface across which energy and momentum are conserved. The equations, given 
in system (|5.35|) - (|5.36|) below, describe the simultaneous time evolution of the 
shock position together with the TOV metric inside the Black Hole such that 
the resulting metric matches a given k = FRW metric Lipschitz continuously 
across the shock, and such that conservation of energy and momentum hold at 
the shock. The conservation constraint, given in 1)5. 25() below, determines the 
TOV pressure, and this is used to close the equations, c.f [H]. Equations (|5.35|) - 
(|5.36|) then determine all other unknowns in the TOV metric beyond the shock, 
and guarantee conservation, once the FRW metric is assigned. In particular, the 
equations guarantee that there are no delta function sources produced by the 
(second order) Einstein equations due to the lower order smoothness (Lipschitz 
continuity) of the metrics at the shock, c.f. (i)-(iv) below. The success of the 
method relies on the fact that once Lipschitz continuity is imposed, the single 
conservation condition (|5.25l) alone guarantees the two conservation constraints 
of the Rankine-Hugoniot jump conditions at the shock, c.f. JHI) and Theorem |21 
below. 

Solutions of equations 1)5. 35(1 - 1)5. 36|) are formally time-reversible without the 
imposition of an entropy condition. In the next section we formulate an entropy 
condition that agrees with the entropy condition of gas dynamics, and we use 
this to construct exact solutions of ()5.35|) - ()5.36|) that describe a class of time- 
irreversible entropy satisfying shocks in which the TOV density and pressure satisfy 
the physical bounds < p < p. 
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To model the expanding universe with a spherical shock wave at FRW position 
f = f{t), we assume that the FRW metric lies in the region f < f{t), and the TOV 
metric lies beyond the shock wave at f > f{t). The corresponding matching for a 
standard TOV metric outside the Black Hole, (and inside the Hubble length), was 
accomplished in jl4 | ll5 | [^. In this section we do not use an EF type regularization 
of the TOV metric, (c.f. H2.9|) ). but rather we employ standard Schwarzschild 
coordinates, avoiding the singularity at j4 = by working either inside or outside 
the Black Hole separately. 

Our proceedure for shock matching is as follows: We first identify the shock 
surface across which a A; = FRW metric matches a TOV metric inside the Black 
Hole, such that the matching is Lipschitz continuous across the interface, and such 
that we have a smooth matching of the spheres of symmetry, (c.f., Lemma 
9, equation (5.3)). The Lipschitz matching of the metrics is guaranteed by a non- 
characteristic condition that always holds outside the Black Hole, when ^ > 0, c.f. 
(|5.2.S)) below. Given this matching, we then determine a conservation constraint 
that guarantees that the Rankine-Hugoniot jump conditions 



hold across the shock. (Here, as usual, [•] denotes the jump in a quantity across 
the shock interface, are the covariant components of the normal vector n to the 
shock surface, and assuming spherical symmetry, we need only require (|5.1|) for 
1^ = 0, 1.) We can then apply Lemma 9, which states that for metrics matched 
Lipschitz continuously across the shock, with a smooth matching of the spheres of 
symmetry, the Rankine-Hugoniot jump conditions (|5.1|) imply that the following 
equivalencies are also valid at the shock: 

(i) The extrinsic curvature is continuous across the shock, 

(ii) The Riemann and Einstein curvature tensors, viewed as second order op- 
erators on the metric components, produce no delta function sources on the shock, 

(iii) In a neighborhood of each point on the shock surface there exist coordinate 
transformations such that in the new coordinates, all second derivatives of the 
metric components are bounded a.e., and 

(iv) At each point on the shock surface there exist coordinate frames that are 
locally Lorentzian. 

In the papers |13| I14| \17\ I21j , the main idea for carrying out this procedure in 
the case ^ > 0, was to show that the single condition 



alone implies the jump conditions ()5.1() . and hence (i)-(iv). However, we must 
modify this idea in the case A < because it turns out that the condition (|5.2j) has 
a non-physical solution which everywhere violates the non-characteristic condition 
used to guarantee the Lipschitz matching of the metrics at the shock. For this 



[T^n^. = 0, 
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reason, to obtain the conservation condition when the shock wave Hes inside the 
Black Hole, we verify the Rankine-Hugoniot conditions directly. 

To start, we match the FRW to the TOV metric when k = and ^4 < 0. The 
A; = FRW metric in the usual (t, r) coordinates is given by, 

ds^ = -dt^ + R{tf [dr^ + r^dn^y (5.3) 

In (|2.18|) - (|2.2fl|) we showed that the mapping (t, r) — > {f,t) that takes metric (|5.3|1 
to standard Schwarzschild coordinates in which the metric takes the form, 

ds'^ = -BpRwdf + A'p.]^y^^df'^ + f^dO^ (5.4) 

is given by 



r = R{t)r, (5.5) 

and 

di = {i)C)dt - {'iljE)dr, (5.6) 

where 

C = (5.7) 

E = -rH. (5.8) 

(Recall that ^ is timelike in the coordinate system (f , t) and spacelike in the 
coordinate system {t,f) when ApRw < 0, H^^ > 1.) Here ip can be taken to be 
any solution of 

^_{i;C) + ^^im = 0, (5.9) 

determined from initial data on any non-characteristic surface, (c.f. the paragraph 
following p.2H) ). Under this coordinate transformation, we obtain 



I _ I 1 

^R'^-f'^R'^ ~ ^ 1 - {rHf 

c.f. ([TTHl) - (fT^ . We use the subscript FRW to distinguish the FRW metric 
coefficients ApRw = 1 — ^^'^^^^ and -B_f_rvk from the TOV metric coefficients 
A = 1 — ^M- and B, that appear in the TOV line element 
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ds'^ = -BdP + A'-^dr'^ + r'^dn'^, (5.12) 

in standard Schwarzschild coordinates. By shock matching, using 1)2. 7|1 and 1)5. lUp . 
we find that 

MFRwit,f) = lf^H^ = ^p{t)r\ (5.13) 
2 D 

Assume now that the (f, t) coordinates that describe the TOV and FRW met- 
rics in standard Schwarzschild coordinates, actually represent a single coordinate 
system for the pair of metrics matched across a shock surface where the metrics 
agree; that is, where A = Afrw and B = BpRw- Setting ApRw = A and using 
(|5.13|) gives 



= (mf = = ^, (5.14) 

r 3 r 

from which we deduce the following conservation of mass condition that must hold 
at the shock surface, and which is independent of V' '■ 

'^p{t)f^ = M{r). (5.15) 

Equation (|5.15|) implicitly defines the shock surface r = f{t), which then determines 
the position r = r{t) = of the shock in the original FRW coordinates {t,r). 
By (|5.14|) we also see that 



N 



I- A 



— = {rHY, 
r 



(5.16) 



and so 



r = y/NH-^ (5.17) 

holds on the shock surface. Using (|5.17|) . we can interpret from the TOV metric 
as follows: if (t, f) is a point on the shock surface, then yjN{t, f) is equal to the 
number of Hubble lengths from the center of the FRW spacetime to the shock 
surface as measured at fixed time t on the FRW side of the shock. Since A > 1 
if and only if ^ < 0, it follows that the shock wave lies beyond one Hubble length 
from the FRW center when ^ < 0. Note that A^ is a convenient variable because 
it appears in system H4.16() - H4.18|) . 

Setting BpRw = B at the shock surface and using (|5.11j) . we obtain the initial 
condition for the integrating factor ip of the FRW metric on the shock surface 
(j5.15|) : namely, 

V; =1 {(1 - i/V)B(f)}"^^^ =1 {AB}-^^^ , (5.18) 
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where we have used that 1 — H^f^ = A < at the shock. Note that AB > 0, and 
the choice of sign in (|5.18)) determines the time orientation for t, c.f. (|5.6|) . 

We now derive a condition that guarantees that the surface f = f(t) is non- 
characteristic when ^ < 0.^'^ To this end, differentiate (|5.15|) with respect to FRW 
time t to obtain, 

M'f = -pr^ + -pf^f. (5.19) 
6 2 

Using (|2.8j) and (|4.17j) in (|5.19j) we obtain the following formula for the speed f, 
the speed at which the shock wave is receding from an observer fixed at the FRW 
origin: 

f = P^Hf. (5.20) 

P + P 

For future reference we record that this directly implies 

Hf f p — 

On the other hand, in light of (|2.21jl . the characteristics for (|5.9|) satisfy 



(5.21) 



f = - {fHY^ + fH. (5.22) 
Thus the shock surface is characteristic if and only if the two speeds in 1)5. 2U() and 



(5.23) 



22|1 are equal, which holds if and only if 

P + P^ 



,p + p} {HrY 

Using ()5.14|) in (|5.23|) we conclude that the shock surface is non-characteristic for 
the PDF if and only if 

A^P±^. (5.24) 
P-P 

The inequality (|5.24j) is not in general ruled out on a timelike shock surface when 
A <Q. However, (|5.24() immediately implies the following lemma which we record 
for use in Section 6: 



^"This is no moot point. Indeed, in previous work, (c.f. ^1 El)) the authors used 
H5.2|l to obtain the conservation constraint across the shock surface when A > 0, but 
we must be careful here because a calculation shows that one of the solutions of (j5.2|l is 
everywhere characteristic for (|5.9|l when A < 0. It follows that H5.1|l does not hold across 
the characteristic surface, and so there is a solution of (|5.2|l that does not represent an 
actual weak solution of G = nT when ^ < 0. It is for this reason that we go directly to 
H5.1(l to construct the conservation constraint in the case A <Q below. 
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Lemma 1 // /), p, p are all positive on the shock surface f = f[t) defined implic- 
itly by i5.15\) . then the surface is non- characteristic for the PDE (1 5. .9)) so long as 
p > p and A < 0. 

Assuming that ()5.24|) holds, the above procedure defines a mapping (t, r) 
(f, t), defined in a neighborhood of any point on the shock surface (f(t),t), such 
that, under this coordinate identification, the FRW metric matches the TOV met- 
ric Lipshitz continuously across the shock surface. Note also that as in the OS 
solution, the coordinate f is timelike in the Schwarzschild coordinates (f, t) when 
^ < 0, but spacelike in the FRW coordinates (t, f), where the coordinate f mea- 
sures arclength distance in each spatial slice t = const. 

It remains to analyze the Rankine-Hugoniot jump conditions 1)5.11) . The fol- 
lowing theorem gives a formulation of the conservation constraint that is amenable 
to analysis. 

Theorem 3 Assume that the coordinate mapping (t, r) — > (f, t) defines a Lipshitz 
continuous matching of an FRW metric to a TOV metric inside the Black Hole, 
in a neighborhood of a point P across a non- characteristic shock surface f = f{t) 
defined by i5.1,5\) . so that that N > 1 at P. Assume further that p, p > at 
P. Then the Rankine-Hugoniot jump relations i5.1\) together with the equivalent 
conditions (i)-(iv) listed after \5.1\) ) all hold at P, if and only if the following 
single conservation constraint holds at the point P : 

^ + \p--p) N 

Note that (|5.25j) immediately implies that if p > p, then also p > p, consistent 
with an outgoing explosion of an inner FRW metric into an outer TOV metric. 
(The construction of an explicit example of such a shock wave is the topic of the 
next section.) Solving for p, it follows that condition ()5.25|) is equivalent to the 
condition 



{p + p)p+{p-p)Np 



(p + p) + {p - p)N 



Note also that from 1)5.26(1 we see that if p = p = 0, then also /O = 0, the conditions 
of the OS solution. 

Before giving the proof of Theorem EJ we first use the conservation constraint 



.26|) to derive the equations that describe the time evolution of the shock surface 
and TOV metric inside the Black Hole that matches a given A; = FRW metric 
at a shock wave positioned beyond one Hubble length > 1. To start, assume 
that p{t),p{t) and R{t) are known functions of the FRW time t that determine a 
unique equation of state p{p), and let a{t) = p{t)/p{t). The matching condition 
()5.15|) written in the form 
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gives p as a function of and r at the shock, and substituting this into the 
conservation constraint (|5.26j) and using p = p{p) gives p as a function of p, N and 
f. The first two equations (flTT^ . (|I?T7)) of the TOV system (|I?TC)l - (|I?TH)) then close 
upon substituting the resulting expression p = p{p, N, r) for p in the first equation 
()4.16|) . The resulting system, 



= i±mm^^. ,5.28, 



N' = -{^ + Kpf\, (5.29) 



then forms a non-autonomous system of two equations in the unknowns (p, N) as a 
function of the independent variable f, the shock position. (Again, "prime" denotes 
^.) Solving (|5.28|) . (|5.29|) gives p{r) and N{r) subject to the initial conditions 



p = po, N = No, at r = f^. (5.30) 

Assuming the FRW solution p{t), p{t), and R{t) is given, we can obtain f as a 
function of t from (|5.27|) . whereby we conclude that p{f) and N{f) determine the 
entire shock wave solution. For example, to set the position of the shock at present 
time in the FRW metric to be at f = fg, choose po = present density, define 



No = "^Porl (5.31) 

(c.f. use H5.27() ). This leaves po as a free parameter. Once we know the solution 
{p{f),N(f)), we can use (|5.27j) to obtain the FRW density /> as a function of f, 
and inverting the relation 



Pit) = 5.32 

gives the shock position f = r{t), the distance from the shock wave to the FRW 
center at FRW time t. The FRW shock position is then 
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The only other unknown in the problem is the metric coefficient B{r) from the 
TOV metric, which we get by integrating the equation (|4.18j) ''^^. 

Bir) = i^oexp {-£ + .m) dc]. (5.33) 

Note that for an outgoing shock wave in an expanding FRW metric, it follows from 
()5.29|) that N decreases and f increases in forward FRW time, and so ^ oo, 
f ^ would correspond to the Big Bang from the FRW point of view. 

We conclude that the problem of constructing FRW-TOV shock waves with 
> 1 reduces to the analysis of the 2x2 system (|5.28p , ()5.29p . The next theorem 
shows that when we make the change of variable p ^ u = ^ and take N to be the 
independent variable, the resulting equivalent system has the nice property that 
the equations are coupled to the FRW metric only through the variable = p- 

Theorem 4 Under the change of variables 

u = V = -, a = -, (5.34) 
P P P 

system i5.28\) . \5.2y\) with conservation constraint i5.2b}) transforms to the equiv- 
alent system 



du r (l + ii) \ \ {2,u - l){a - u)N + Qu{l + u) 



dN l2(l + 3n)iVj I {a-u)N + {l + u) 
df If 



with conservation constraint 



dN l + 3uN' 



-a{l + u) + {a - u)N 
(l + u) + {a -u)N 



(5.35) 



(5.36) 



(5.37) 



For such solutions, the speed of the shock interface relative to the fluid comoving 
on the FRW side of the shock, is given by 

s = Rr = ^ (tT^) ■ ^^'^^^ 



"'^^Note that (|5.33(l completes the definition of the TOV metric, and so the integrating fac- 
tor -0 for the shock matching is determined by H5.18|l . and from this the non-characteristic 
condition (|5.24l) imphes the Lipschitz matching of the metrics. 
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Note that the dependence of (|5.35() - H5.37|) on the FRW metric is only through 
the variable a, and hence (|5.35j) is coupled to (|5.36j) only through the function a. 
In particular, if we assume that a is constant in the FRW metric, (a reasonable 
model problem, c.f. |14j). then equation (|5.35|) uncouples from equation (|5.36|) to 
form a non-autonomous scalar equation in u and which is amenable to phase 
plane analysis. This is the starting point for the exact solutions constructed in the 
next section. 

Note too that since (OHl) implies that ^ < 0, and VN gives the number 
of Hubble lengths to the shock wave, c.f. (|5.17jl . equation (|5.3H|) implies that the 
number of Hubble lengths from the FRW center to the shock wave decreases in 
time, as claimed in the introduction. 

Since system system 1)5. 35() . ()5.36() is equivalent to system H5.28() . H5.29() . the 
former system also determines all quantities in the shock wave solution, but now 
as a function of the variable A^. For example, a solution of system (|5.35|) . ()5.36|) 
is determined by the initial conditions u = uq and f = fo at N = Nq, which 
entails assigning three constants. Alternatively, if we ask that the shock wave be 
positioned \/Nq Hubble lengths from the FRW center at the present FRW value 
of the density p = po, then by (|5.27jl , the initial shock position is determined by 



ro = W , (5.39) 

V i^Po 

and there is then a one parameter family of such solutions {u{N),f{N)) determined 
by the remaining free parameter 



no = ^. (5.40) 
Po 

Once we know the solution (u{N),f{N)), we can use 

m 



K [r{N)] 



2 ' 



(5.41) 



to obtain p{N), the FRW density p as a function of A^ at the shock. Since the 
FRW metric is known, knowing p{N) determines p{N) and o"(A^) = , and 
using these in (|5.37)) we obtain v{N) and p{N) = v{N)p(N). Finally, the TOV 
metric coefficient B is determined as a function of A^ by integrating the equation 
(|4.18)1 with respect to A^, yielding the formula 



B(,) = B„ exp I - £ ^ (^ + .p«)) A . (5.42) 



32 



To connect N with the FRW time t, invert the relation H5.41() . obtaining N{t). 
Since R(t) is assumed known, we can use 

to get the FRW shock position r as a function of t. Relations (|5.39|) - (|5.4.Sj) formally 
determine an exact shock wave solution of the Einstein equations for any given 
{k = 0) FRW metric, and any given solution of system (|5..S5j) . (|5..S6j) . 

Proof of Theorem |4j Equation (|5..S7j) follows directly from (|5.26j) upon dividing 
through by p and making the substitutions in ()5.34() . Equation ()5.36() follows 
directly from (|4.17|) upon noting that Kpf = a consequence of (j5.32|) . From 
we have 



1 _ IN' 
f ~ ^(l + 3n)77' 

To verify (|5.;-i5|) . start with 



Then by (1^:^ and (t07|) . 

p' u + v N' If {l + u){a-u) 



p 2 N-1 2 {{1 + u) + {a -u)N 
and using (|5.27|) . 



Kpr"^ J r J \ 1 + 3u J N ' 



where we have used (15.441). Thus 



and so using (|5.4H|) and (|5.48|) in (|5.45j) gives 



u' 



If {l + u){a-u) \ ^, f3u{l + u)\N' 



2 {{l + u) + {a- u)N J \ l + 3u J N 

(a - u){3u'^ + 2u- 1)N + 2u{l + u)'^ 
2(1 + 3u)[{l + u) + {a- u)N]N ' 
{l + u)N' \ f {3u-l){a -u)N + 6u{l+u) 



2(1 + 3u)N J i {a- u)N + (1 + n) 



(5.44) 



pj p p^ 



N', (5.46) 



p , ( 3u{l + u) \ N' 

7^ = l^T3^j^' 
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which upon dividing by N' verifies ()5.35|1 . Finally, since the fluid on the FRW side 
of the shock is assumed to be co-moving with the radial coordinate r of the FRW 
coordinate system, it follows that s = Rr gives the speed of the shock relative to 
the FRW fluid, (c.f., the discussion after Theorem 4 of and so (|5.38|) follows 
directly from (|5.21|) . □ 

Proof of Theorem |3J Assume that we are given an FRW metric and TOV 
metric inside the Black Hole that match Lipshitz continuously across a smooth, 
non-characteristic shock surface f = r{t) defined in a neighborhood of point P, 
such that the hypotheses of Theorem |31 hold. We show that ^b.\^ holds if and only 
if (|5.25|) holds. We use the following lemma: 

Lemma 2 On the shock surface, in (f, t) coordinates, we have. 



[T] 



\rp^U rpflV 1 

i^FRW ~ ^TOV\ 

(p + p)N + {p + p){l 
ipVN{p + p) 



(5.49) 



and 



no = ^p(N^-^ + l), (5.50) 
\ p + p J 

where n^dx^ = fiodf + nidi is the covariant normal to the shock surface. 

We use Lemma 121 to verify Theorem El and postpone the proof of Lemma |2l until 
the end. To verify that (jS-lj) is equivalent to ()5.25() at a point P on the shock 
surface, it suffices to show that det{[rY'^) = holds at P if and only if (|5.25j) 
holds at P, and that when this holds, the shock normal n is in the kernel of [T]. 
Indeed, if ()5.25|) holds at P, and we know n is in the kernel of [T]^'', then 1)5. 1|) 
and (|5.2|) hold at P, and so our general theory in [T^ would imply that (i)-(iv) 
hold at P as well. Conversely, if (i)-(iv) and ()5.1|) all hold at P, then the RH jump 
conditions (|5.1|) alone imply that det{\TY^) = holds at P, and hence we would 
have ()5.25|1 at P as well. 

We first show that det ([T]^'") = holds at P if and only if holds at P. 

But by dsani), 



det mn = {{p + P)N + {-p + p){l- N)}, {ilj\p + p) + {p- p)N]] 

-ip'^N{p+pf = {), (5.52) 

and solving for p gives 
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p 



N{p + pf 



{{p-p)N + p + p}j 
{{p + pN){.},-N{p + pf] 



+ {p + pN) 



N -I 



(iV-l 



{N- 
{p{p- 



-p)N- 



p{p + p)} 



{N-l){{-p + p) + {p-p)N} ' 

from which H5.25() fohows at once. (Note that we have assumed without loss of 
generality that {•}/ 7^ 0, which is valid because, if {-^i = 0, then (|5.52j) implies 
p = p = 0, and /9 = violates N > I, e.g., ^ 7^ 0, c.f. (f^3^ . (EUl).) Thus it 
remains only to show that if det {\TY'^) = holds at P, then the shock normal n 
at P is in the kernel of [T]. To this end, note that by 1)5. 5U() and 1)5. 51(1 . 



N 



P-P 
P + P 
P + P 
P + P 



+ 1 ) ^/N{p + p) 



(5.53) 



{{p + p) + {p-p)N 



and so n^\TY = as well because we assume (iet[T] = 0, and [T] is a 2 x 2 matrix 
The proof of Theorem |21 is complete, once we give the 

Proof of Lemma |2j To verify (|5.49j) . we calculate Tj^qv '^frw^ stress 
tensors in barred coordinates x = (f, t), at the shock. (Recall that x'' = f is the 
timelike coordinate, and = i is spacelike because > 1.) Using the formula 



Ttov = (P + p)utov^tov +P9'"', (5.54) 
and the assumption that the TOV fluid is co-moving in barred coordinates, 







VN -1 











we obtain 



-'-TOV 



(p + p) 



-A 









+ pA 









^2 



^1/ 



-'^^ Although H5.53|l appears to be a miracle here, this is to be expected from where 
in the case A > 0, the single condition [T'^'^Jn^rii, — alone implied the equivalencies (i)- 
(iv), as well as the Rankine-Hugoniot conditions (|5.1|l . One also expects this in light 
of the fact that the Einstein equations for a spherically symmetric metric in standard 
Schwarzschild coordinates, contains only one second order, (i.e. for Lipschitz metrics, one 
weak), equation, so we would expect only one jump condition. 
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Similarly, 



'^FRW = iP + P)uFRW^FRW + P9 



111/ 



(5.55) 



but in this case the assumption that the FRW fluid is co-moving in FRW unbarred 
coordinates x"' = [t,r)'^, implies that 



where 



u 



dx"' 



1 




(5.56) 



(5.57) 



gives the components of the FRW fluid velocity in standard FRW coordinates 
(t,r). 

Claim: When A; = 0, we have 



dx^ 



Vn r 



(5.58) 



Proof of Claim: The coordinate mapping r = Rr implies 



= Rr = Hr, 

at 



dxP 
dr 



R, 



which verifies the the first row of ()5.58|) upon using H5.16() . Hr 
by 



Substituting 



into (|5.60j) gives 



df = d{Rr) = Rrdt + Rdr = Hfdt + Rdr 



(5.59) 
N. Moreover, 

(5.60) 



= V {dt + HfRdr) 



(5.61) 
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from which we can deduce the second row of H5.58() . which verifies the Claim. 
Continuing now from (|5.56jl : 



^ R 



- M 


' 1 ' 


a 




- a 










Using this in H5.55() gives 



TpRW = {P + P) 
Using (|5.H2|) together with (|5.55|) gives 



N TpVN 



1 

0-^2 



{p + p)N + {p + p){l - N) i^^ip + p) 

^I;VN{p + p) i:^{{p + p) + {p-p)N} 



which is (12311). 

Finally, we use the Claim to verify 1)5. 50() and ()5.5ip as follows. Let 



(5.62) 



(p{t,r) =r -r{t) = 0, (5.63) 
where r(t) denotes the shock surface in FRW {t, r)-coordinates. Then 

dcj) = —rdt + dr = n^dt + nidr, 
gives a covariant normal nadx" to the shock surface in x = (i, r) coordinates, 



no 
ni 



1. 



(5.64) 



In X = (r, f)-coordinates, 



(5.65) 



where ffjr is the inverse of the matrix (|5.58|) . 



dx" 



dxi" det 



dx°' 



tPRVN - R 



(5.66) 
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Using ()5.64|) and ()5.66|) in ()5.65|) and neglecting the irrelevant factor , , , we 
obtain, (to within a minus sign), 

= ^ \rR^ + l] , -rR - ^iV ) , (5.67) 
and substituting H5.21() in the form 



N p — p 
it p + p 

into (|5.67|) . then gives the formulas for given in (|5.5())) . (|5.5H) . The proof of 
Lemma 121 is complete. □ 

We summarize the results in this section as follows: System (|5.35I) - (I5.36() de- 
scribes the TOV metrics that match a given k = Q FRW metric across a shock 
wave discontinuity when A < 0, assuming that the conservation constraint ()5.37|) is 
satisfied. Solutions of these equations, together with the initial conditions in (|5.:-{9j) - 
()5.43|) . determine weak solutions of the Einstein equations containing a shock wave 
interface across which the metric is only Lipschitz continuous, and such that con- 
ditions (i)-(iv) (after 1)5. 2|) ) are satisfied. Since the resulting solutions are formally 
time-reversible, it still remains to impose a physically meaningful entropy con- 
dition. The entropy condition for shocks determines the time orientation of the 
solution, c.f. ^2j. In the next section we construct a class of exact solutions of 
these equations which satisfy an entropy condition that agrees with the entropy 
condition of gas dynamics in the non-relativistic limit. 

6 Exact Shock Wave Solutions When A < 0. 

In this section we derive a class of exact solutions of equations ()5.35l) - ()5.36() in the 
case when the FRW pressure is given by the equation of state 



p = ap, (6.1) 
where a is assumed to be constant. 



< cr < 1. (6.2) 

In this section, as an entropy condition, we impose the condition that the shock be 
compressive, a condition sufficient to choose the physically relevant stable shocks 
in classical gas dynamics, jU [71 ^J. That is, we will require that the pressure 
and density be larger on the side of the shock that receives the mass flux. Exact 
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solutions satisfying the equation of state ()6.1() were constructed in for the case 
A > 0, so we can interpret the results here as an extension of the results in ^Ij 
to the case A < 0. However, in the conservation constraint led to a TOV 
equation of state also of the form (|6.1|) . but in our case here, the TOV pressure is 
not so simple due to the fact that the conservation constraint (|5.l-{7|) is of a different 
form. 

Assuming the solution to the FRW equations (|2.7j) . (|2.8|) is given in 

Theorem 121 equations ()2.41|) - (|2.45|) . and we assume an expanding universe, {R > 
0), with initial conditions R = at t = 0, and R = 1, at t = to. The solution has 
one assignable constant which can be taken to be to, po or Hq, since 

4 1 . , 

Po = o..r, I (6.3) 



Ha 



3^(1+^)2*2' 



3(1+ a) 



^, (6.4) 
to 

where we interpret to ^-s "present time" in the FRW metric^^, and the relations 
(|09ll - (l5:i3]) hold. 

We now solve the TOV shock equations (|5.35j) - (|5.3(ij) . to obtain a shock wave 
solution that matches this FRW metric for all R{t) in the interval between R = 
and R = 1. Since (|5.35|) is coupled to (|5.36j) only through the function a, our 
assumption that a is constant in the FRW metric implies that solutions of system 
(|5.35j) - (|5.3(i|) are determined by solutions of the scalar non-autonomous equation 
()5.35|) . Making the change of variable 5 = 1 /N, (we do this to transform the "Big 
Bang" N ^ oo over to rest point at S" ^ 0, c.f. [25), equation (|5.35|) takes the 
form 

du _ f (1 + n) \ f {3u-l){a-u)+6u{l + u)S \ 

(i5 ~ l2(l + 3n)5j I (cT - n) + (1 + n)5 J' ^ ' 

where u = p/p. We now construct solutions of (|6.5() that model the "Big Bang" as 
a localized explosion with an outgoing blast wave emanating from f = at time 

^^Since we are assuniing the idealized equation of state H6.1|l for the FRW metric, the 
value of the Hubble constant at present time alone determines the solution. In the standard 
model of cosmology based on a A; = FRW metric, Pj, the solution is determined by two 
assignable constants, the Hubble constant Hq and the background radiation temperature 
To, Ql^. Assuming (|6.1I) . the freedom to assign Tq is represented by the freedom to 
assign a. Also, in the standard model the galaxies follow the particle paths of the matter 
field after the time when the pressure of the matter field is essentially zero. In the model 
H6.1|l . the matter is co- moving with respect to the radial coordinate r, and hence the 
particle paths follow r = 0. Loosely speaking, we refer to the motion along the particle 
paths as the motion of the galaxies in the cosmological interpretation of a general FRW 
metric. Thus we say that the shock wave is exploding outward through the galaxies when 
the shock wave satisfies f > 0, 
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t = 0. Thus, motivated by classical gas dynamics, we seek solutions of (|6.5|) that 
satisfy the entropy conditions, 



0<p<p, (6.6) 
0<p<p, (6.7) 

and meet the physical bounds on the TOV equation of state 

0<p< p. (6.8) 

Conditions 1)6. 6|) , 1)6. 7|) for outgoing shock waves imply that the shock wave is com- 
pressive. The condition ()6.8j) implies that the TOV equation of state is physically 
reasonable. Note that the conditions > 1 and < p < p restrict the domain 
of ()6.5|) to the region < u < a < 1,0<S'< 1. The following theorem shows 
that the inequalities 1)6. 6(1 . (|6.7() and 1)6. 8|) can all be expressed in terms of u and 
S alone, and are all implied by a single inequality. 

Theorem 5 Assume that 

0<u<a, (6.9) 

< 5 < 1, (6.10) 

and the conservation constraint (|5.37| ) holds. Then the bounds \6. 6\) - W^) are all 
implied by the single condition 

S<C^)(^^)^E(u). (6,11) 



1 + uj \a + u 



Proof of Theorem |5} Multiplying (|6.9)) through by p > gives ()6.6)) . To verify 
that p < p, write 1)5. 37() in the form 



S + a ' 
where 

a — u 



^^ + " (6.12) 



1 + u 



> 0, (6.13) 



and observe that p < p is equivalent to f < 1, and by 1)6.12(1 . v < 1 ii and only if 
—a < 1. To verify ()6.8|) . observe that p < p is equivalent to ^ > 1, and by ()6.12j) . 
this is equivalent to 

—aS + a ^ 
{S + a)u ' 
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which is easily seen to be equivalent to ()6.11|) . Thus it remains only to verify < p. 
But observe that p > is equivalent to t> > 0, and by (|6.12j) . this is equivalent to 

S < (6.14) 

a 1 + u 

which is implied by (|6.11j) because 

E{u) < 



a 1 + u 

The proof of the theorem is complete. □ 

6.1 Analysis of the equations 

The purpose of this section is to prove the following theorem: 

Theorem 6 For every a, {) < a < 1, there exists a unique solution Ua{S) of I ^.,5)) . 

such that / 1 6'. y\) and i6.11\) hold on the solution for all S, < S < 1, and on this 
solution, 

< u^{S) < u, (6.15) 

lims^o Ua{S) = u, (6.16) 

where 

u = Min{a, 1/3} , (6.17) 

and 

lim^^ip = = lims_+i p. (6.18) 

Proof of Theorem |5J Solutions of 1)6.5(1 are determined by trajectories of the 
autonomous system 

S' = 2S{l + 3u){{a-u) + {l + u)S} = F{S,u), (6.19) 
u' = il + u){-il-3u)ia -u) + 6uil + u)S} = G{S,u), 

(6.20) 

which has two degenerate rest points, S = 0, u = a and S = 0, u = 1/3, in 
the (5*, tt)-phase plane. Here "prime" denotes differentiation with respect to a pa- 
rameter ^, and we recover equation (j6.5j) by eliminating ^ via du/dS = u'/S'. 
Since -E(O) = 1, (c.f. H6.11() l. the "initial condition" no-(l) = is a consequence 
of the entropy inequality (|6.11j) . and thus the uniqueness of u„{S) follows from 
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uniqueness of solutions for system H6.19() . H6.2U() because u = 0, 5* = 1, is a reg- 
ular initial condition for the system. Thus (|6.18j) follows from (|5.34|) and (|5.37|) 
because lim5^i Ua-{S) = 0. It remains to prove the existence of the solution Ua{S) 
on < 5 < 1, and to establish (|6.16|) . 

Our analysis of these equations is based on the construction of invariant re- 
gions. Before we analyze the phase plane for system ()6.19|) . ()6.2U|) . first note that 
F{S, u) > in the region 0<n<(T, 0<5'<1, and the isocline G = is given by 

{a -u){l/3-u) _ 
^= 2u{l + u) ^'^^^'^^^ ('-'^^ 

Let Ra denote the region^^ 

i?^ = {(S, n) : < n < n, < 5 < h{u, a)} , (6.22) 

c.f. Figure 1. 



Lemma 3 The region R^j is a negatively invariant region for system \d.iy\) . \6.2U\) 
for each a between and 1. 



Proof: To show that R^ is negatively invariant, it suffices to show that the vector 
field (F,G) points out of or tangent to R^ on the boundary of R^, as shown in 
Figure 1. Now the vector field (F, G) is tangent to the axis S = because S* = is 
a solution trajectory of system 1)6. 19() . ()6.2U|) . Since the boundary u = lies below 
the G-isocline, it follows that G < along the boundary u = 0, and so the vector 
field {F, G) points out of R^ along the lower boundary. Finally, the vector field 
points out of Ra along the upper boundary S = h because, for each a, < a < 1, 
S = h is a G-isocline, F is positive, and S = h{u, a) is a decreasing function of 
u on < u < u. Indeed, S = h{u, a) is decreasing on < u < u because, (using 
()6.21|) . and prime for 



where ^ = ^ and 

A = (o--n)(l/3-u), 
B = 2n(l+n), 

and we use the fact that A > 0, A' < 0, B > 0, and B' > in < u < u. 

^^We include the line m = in Ra- in order that the initial point u — 0, S — 1 for the 
solution trajectory lies in R„. However, the entropy conditions H6.6() . (|6.8() do not hold at 
the transitional point u = 0, S" = 0, the place where the solution Ua-{S) continues naturally 
to the OS solution when p << 1. 
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Now the existence of the orbit Ua{S) on the entire interval < S* < 1 is a 
consequence of the fact that S = 0, u = u is the only rest point in the region R^^. 
That is, the unique trajectory of system (|6.19() . (|6.2U() starting from initial point 
{S,u) = (1,0) G Ra, must tend to the unique rest point (0,u) in backward time 
^ — > — oo. Moreover, since G < and F > in the interior of Ra it follows that S 
and u are monotone along the orbit, so the orbit defines the trajectory u = Uo-(S'), 
< 5 < 1, as well as its inverse, S = Sa{u), < u < u. 

It remains only to verify that the entropy condition H6.11() holds all along the 
solution, and for this it suffices to show that 

A(n) = E{u) - Sa{u) > 0, (6.24) 

for all u£ [0,u). But 

A(0) = (6.25) 

because E{0) = 1 = 5(j(0), and 

A{u) > (6.26) 

because E{u) > and Sa{u) = 0. Thus to verify that A(n) > for < u < n, it 
suffices only to show that A{u) > in a deleted neighborhood of u = 0, and that 
A'(u) = has at most one root in < u < u. But differentiating (|6.11j) gives 



E'(u) = -2(1 + a) 



(T — 



((j + 'u)2(l+n)2- 
Using this and (|6.5() at u = gives 



and 



E = 4 2 — > -6 2 — = '^'^^)' 



which implies that A'(0) = 0, A"(0) > 0, so A(u) > 0, near u = 0. Moreover, by 
(|6.1H) and (|6.5|) we also see that A'(ti) = if and only if 



52{(a + n)2(l + u)3(l + 3n)}^ (6.27) 
+S{{a- u){a + uf{l + u)2(l + 3n) + 6n(l + cr){l + u)'^{a - n^)}^ 
+ {(3n - 1)(1 + a){a - u){l + u){a - u^)} = 0. 

But j-}^ > 0, j-}^ > 0, and j-}^ < for < u < u because u = Min{a, 1/3} , and 
thus the quadratic (|6.27j) has exactly one positive root. We conclude that (|6.24|) 
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holds. Finally, to verify H6.18() . we note that u ^ as S" — > 1, so p — > 0. The 
fact that liuis^ip = follows directly from (|5.,S7|1 . The proof of Theorem 1^1 is 
complete. □ 

Theorem ini implies that the entropy conditions ()6.6() . ()6.7() pick out the unique 
solution of system (|6.5|) emanating from the initial point S = l,u = ^ = 0= p. 

6.2 The shock speed 

In this section we estimate the shock speed Sa{S) = s{ua{S)), the speed of the 
shock relative to an observer fixed relative to the FRW fluid, along the solutions 
Ua{S) of TheoremEl < 5 < 1, < u^{S) < a. 

Theorem 7 Let < o" < 1. Then the shock wave is everywhere suhluminous, that 
is, 

K{S)\ < 1, (6.28) 

for all < S < 1, if and only if 

a < 1/3. (6.29) 

Proof: Formula (|5.H8j) gives 

sAS) = iir = (^) , (6^30) 

and thus since lims^oUa{S) = u = Min{a, 1/3}, it follows at once that lims^oSf^{S) = 
oo if cj > 1/3. Since lims^oUa{S) = a when a < 1/3, we conclude from ()6.3U|1 
that the precise value of the shock speed Sa in the limit 5^0 depends on the 
asymptotic behavior of the solution Ua-{S) as 5 — > 0. This is addressed in the next 
subsection. 

To verify H6.28() . ()6.29() . we first prove the following lemma. 
Lemma 4 Assume that < a < 1/3, and let S = Qa{u) = Q{u) he defined by 

^?aM=«'|^^- (6.31) 

Then the region 

Q''^^{{S,u):S>Qa{u)}nK, (6.32) 

is a negatively invariant region for solutions of system 116. iy\) . \6.2U\) . so long as 
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Moreover, for a = 1/3 and sufficiently small e > 0, there exists a 6 > 0, such that 
if 



a' < (6.34) 
then solutions of i6.iy\) . h6.2U\) . starting in the region 

Qv3(e) = 2V3 n {(5, u) : < < 1/3 - e} , (6.35) 
can only leave the region Q'^i^{e) through the boundary u = 1/3 — e. 

Note first that, assuming a < 1/3, the condition a"^ < guarantees that Qa{u) < 
h{u) in the region 0<u<(T, 0<5<1 because by (|6.21|1 . Qa{u) < h{u) is 



equivalent to 



and since 



< C-I^) <l±iil, (6.36) 

\ a — u J 2u 



2a \ a-u J 2u ' ^ ' 



a 



< 2^ guarantees ()6.36l) in < u < cr, < 5" < 1. It fohows that when < 
the region is the region in the {S, n)-plane that hes above the curve S = Qa{u), 
above the curve u = 0, and below the isocline S = h{u). Moreover, when < 
the initial point u = Q, S = 1 for solution trajectory u^(S) lies within the region 
Q'^ because Qa{0) = a^o'^ < 1. 

Before proving the lemma, we first use Lemma 0] to complete the proof of 
Theorem |7| by verifying ()6.28|) . Since the initial condition S = 1, u = lies within 
the region Q", we can estimate the shock speed all along the orbit U(j{S) by 

1 / a — u\ 1 / 1 + u\ / a — u\ 1 



and so taking a = < (optimal for (|6.33j) ). we obtain 



s„{u)<^ <\, (6.39) 

for all < cr < 1/3, < u < a. This verifies K!M for < < 1/3. However, 
(|6.39|) does not provide a strict inequality if o" = 1/3. 
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To obtain the strict inequality (|6.28|) in the case a = 1/3, we use that the 
initial condition = 1, u = also lies in Q^y3(e) for a = , so if we assume 
Ml/3 (5) ^ 1/3 — e, then (|6..S4j) . (|6..S5j) imply that the shock speed in the region 
Q"/3(e) can be estimated by 



si/3{u) 



< 



VT^ < 1, 



(6.40) 



when fj = 1/3. Since e is arbitrary, we conclude that H6.28|) holds for all < a < 
1/3, < S* < 1. 

The proof of Theorem [7| is complete, once we give the 

Proof of Lemma Ht Assume 1)6. 29() . To show that QJ^ is negatively invariant, 
it suffices to prove that the vector field (F,G) restricted to the curve S = Qaiu), 
always points into the region below the curve. That is, it suffices to prove that 



IQ'(n) 
is 



> 



on = Qa{u) = Q, where 4^ is given in ()6.5() . Thus our condition 



ig'l = o^2 (l+^)(^-^) > 2(1 + 3u)Q[{a - u) + {I + u)Q] 



(l + n)3 - (i + u)\{3u-l){a-u)+6u{l + u)Q\' 
After some algebra, this reduces to 



(1 + 0-) ^ {l + u) + {a-u)a 



2 



{a - n)(l + 3n) " |(1 - 3'u)(l + u) - 6u{a - u)a'^\ ' 

(6.41) 

But the term inside the absolute value is positive so long as a < because then 



< — < 
- 2a - 



1 — 3n\ 1 + u 



a 



u 



(6.42) 



as a result of (|6.37j) . 

Therefore, using that the term in the absolute value in (|6.4H) is positive, we 
can solve 1)6. 41() for to obtain the following condition equivalent to the condition 
that S > Q{u) is negatively invariant: 



< 



(1 + u) [(1 + a)(l - 3m) - (1 + 3u){a - u)] 
(1 + 3m)(c7 - m)2 + (1 + o-)6m(o- - n) 

(1 + 3n) 



fl + n) 



a-3«) 



(1 + 3M)(f7-n) + 6(1 + o-)u 



(6.43) 



But using that 



1 - 3n 



a 



u 



>3, 
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it follows that the inequality 

(1 + 6u)[a — n) + 6(1 + a)u 

guarantees (|6.43j) all along the curve 5 = Qa{u). But (j)'„{u) <OforO<n<(j< 
1/3 because (paiu) = 4 where 

A = {I + u)[2 + 3{a - u)], 

B = {1 + 3n)(cj - n) + 6(1 + a)u, 

and 

BA' - AB' < 0. 

Therefore, 

(t>a{u) > 4>a{(^) = — • 

3a 

From this we conclude that if (|6.33|) holds, that is, if < < then (|6.44j) 
holds, and hence is an invariant region for all < a < 1/3, as claimed. It 
remains to verify (|6.34|) and (|6.35|) in the case a = 1/3. 
So assume a = 1/3. Then from ()6.44l) . we have 

(l + n)[2 + 3(1/3-^)] 
^^^^^ ' (l + 3n)(l/3-u) + 6(l + l/3)7x' 

But since, (j)'^^^{u) < for < u < 1/3, it follows that for every e > there exists 
a (5 > such that 

01/3 W > 

in the region 0<n<l/3 — e. Thus, if we choose < then 1)6. 44() holds, and 
so S > Qi/siu) is an invariant region for < w < 1/3 — e. It follows that (|6.34|) 
implies that orbits can only leave the region Qy^ie) through the line u = 1/3 — e, 
as claimed. The proof of Lemma 0] is complete. □ 

6.3 The shock speed at the Big Bang 

In this section we calculate the shock speed Sa{S) along the orbit u = Ua{S) in 
the limit S" — > 0, the instant of the Big Bang. In Theorem [3 we showed that 
lims^QSfj = oo if o" > 1/3, and the following theorem asserts that lims^oSa = 
for cj < 1/3, but that lims^oSa = 1 for the special value o" = 1/3. In particular, 
this confirms that the estimates of Theorem |3 are sharp. 

Theorem 8 Let < o" < 1/3, and let s^j^S) = Sa{ua-{S)) denote the shock speed 
along the solution u = Ua{S), given in i6.S(jl\) . Then if a > 1/3, 

lim Sa(S) = oo, (6.45) 
5^0 
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if a <!/?,, 



and if a = 1/3, 



lim s^(S) = 0, (6.46) 
5^0 



lim Sa(S) = 1. (6.47) 
5^0 



Proof: By (|6.,3fl|) . the shock speed is given by 



sAS) = -^J^^-^), (6.48) 



from which H6.45() is evident. To verify (|6.46() . we show that when < cr < 1/3, 
and a = there exists m > such that 

V OCT 

(1, 0)GR^^ {{S, u):S>m{a- u)} H Q«, 

and that the region is negatively invariant for solutions of system (|6.19j) . c.f. 
(|6.32|) . If is negatively invariant and contains the initial point {S,u) = (1,0), 
then Ucr{S) > a — 5 in a neighborhood oi u = a, which implies ()6.46|) because 
then we can use (|6.48|) to conclude 

< lim s^{S) < lim (r-^^TTK) = 0- 
To see that is negatively invariant for some m > 0, assume first that 

c.f. (|6.21|) . This implies that ■m{a — u) < h{u) in the region 0<n<(T, 0<S'<1. 
Thus when m < h'{a), the region consists of the set of all points in Q° that 
lie above the curve S = m{a — u) and below the curve S = h{u). It is also easily 
verified that the initial condition S = \, u = Q lies in BJ^ when m < h'{a). 

To find a value of m < h'{a) for which is negatively invariant, we mimic 
the proof of Lemma 0] using L = L{u) = m{(T — u) in place of Qa{u). That is, 
is negatively invariant if the vector field {F, G) points into the region below 
S = L{u), all along S = L{u), which is to say that |i^'(u)| > ^ on = L{u), 

where ^ is given in (|6.5|) . Since we have already verified that is negatively 
invariant for a = —k=i we actually need only verify the negative invariance of the 
curve S = L{u) for values u < u < a, where u is the value of u at the point of 
intersection {S,u) of the line S = L{u) and the quadratic S = Qa{u) = 3^ (1+^)2 ; 
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c.f. H6.31() . Since L is linear and Qa is quadratic at (0, a), it follows that n — > a as 
m ^ 0. Thus our condition is 

' I (I u)\{a - u){'iu - I) + Qu{l + u)Ly 

for u < u < a. After some algebra, this reduces to 

- |(l+n)|(3n-l)+6n(l+n)m| ^ ^ 

But the term inside the absolute value is negative because 

^^^TTTT^' (6.51) 
6u[l + u) 

a consequence of the fact that we are assuming m < h'{a). Therefore, we can solve 
(|6.5U() for m to obtain the following condition equivalent to the condition that 
5 > L{u) is negatively invariant: 

m<- — — , '--=i;a(u). 6.52 

- {l + u)[2{l + 3u){a -u) + 6u] ^ ' 

Now assuming cr < | is fixed, it follows that ^|) = A/B, where 

B = {1 + n)[2(l + 3u){a - u) + 6u] < 2(1 + a)i7a + Sa"^), (6.53) 

and 

A = 1 - 3^ - 2(1 + 3u){a - u). 

Therefore, to verify (|6.52)1 for u < u < a for some m > 0, it suffices only to show 
that for m sufficiently small, A = 1 — 3u — 2(1 + 3u){a — u) is bounded uniformly 
away from zero for all u E [u, a]. But 1 — 3(T = e > 0, so limm^o u = cr implies that 
there exists an (5 > such that m < 6 guarantees that 2(1 + 3u){a — n) < e/2 for 
all u G [u, a]. In this case, A = l — 3u — 2(1 + 3u){a — u) > |, and hence by (|6.53j) . 

It follows from ()6.52|) that if we choose m = Min{6, 6i}, Then for fixed o" < ^, it 
follows that 

m < tpa{u), (6.55) 

for all u £ [u,a], and hence that -R™ is negatively invariant. Thus the proof of 
()6.46|) is complete, and it remains only to consider the case a = 1/3. 

In the case a = 1/3, is not negatively invariant for any m > 0, and to 
evaluate lim5_>o Sfj(S) in this case we need the asymptotics of the solution Ui/^{S) 
near 5 = 0. This is given in the following lemma: 
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Lemma 5 Assume a = 1/3. Then 

UaiS) ~ 1/3 — m^:^/S, as S ^ 0, (6.56) 

where 

4 

= -. (6.57) 

Lemma ini implies (j6.47|) of Theorem |H1 (the case o" = 1/3), because using (|6.56|) in 
(|6.3U|) we obtain, 



4 



^/S\l + Ua{S)J 1 + i 

Note that (|6.56|) confirms that the orbit ui/^{S) comes into the rest point S = 0, 

u = 1/3, asymptotically like the negatively invariant curve 5 = Qi{u) = (| — 
uf, c.f. (lOTTl . 

Proof of Lemma |SJ Assume a = 1/3. To verify (|6.56|) . write 

u^{S) = 1/3 + (piS), 0<5<1. (6.59) 

We find the equation that cp satisfies asymptotically as 5 ^ 0. (Note that (j){0) = 
because Uo-{S) tends to the rest point (0, 1/3) as 5 — > 0.) Putting (|6.59j) into (|6.5|1 
we obtain 



Scl>'iS) 





-3(/)2 + (2 + 6,/.)(| + (/>)5 


2(2 + 3(/)) 







(6.60) 



But by Lemma ^ we know that Q^/^ is negatively invariant for a = 1, and from 
this it follows that ^1/3(6') is squeezed below the G- isocline and above the curve 
u = Q-\S) = 1/3 - VS as S ^ 0. It follows that (1/3 - 7x1/3(5)) is order S^/^ as 
S ^ 0. Using this in (|6.6()|) gives 



50'(5) 



(I) -3<t>' + IS 



2(2) 



(6.61) 



where ~ means to leading order in 5 as S" 
equation (|6.6H) takes the asymptotic form 



0. Thus, to leading order in S, 

(6.62) 
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which we can write in the form 



(0^)' = ^ - (6-63) 
which is hnear in v = (j)'^ with 

a = 2, (6.64) 

and 

P = 16/9. (6.65) 
But (|6.6l'{j) has the general solution 

v = kS'' ^S, (6.66) 

1 — a 

where k is a constant. Since a > 1, the first term in ()6.66|) is higher order in S, 
and thus we have shown that to leading order in S, equation (|6.62)) has the unique 
solution 



/,2 



4>{S)=mWS, (6.67) 

where 



= W^— = 4/3. (6.68) 
V a — 1 

This agrees with (|6.56|) and (|6.57|) . Thus the proof of Lemma 0] will be complete 
once we prove the following lemma. 

Lemma 6 Assume a = 1/3. Then every solution u{S) of 116. 5\) that enters the rest 
point 5' = 0, M = l/3 from inside the invariant region Qi^^, a = 1, must satisfy 

w(S) ~ 1/3- asS^O. (6.69) 

Proof: By the derivation of (|6.62j) . we know that 

n(5) = 1/3 + 0(5), 

where (j) satisfies 

{^^y = a^-P-e{S), (6.70) 



52 



for some function, e{S) — > as 5 ^ 0. Thus it suffices to sliow that hms_»o 
0. But ()6.7U|) is hnear and has the general solution 

02(5) = KS'' + 5" {f3 + e(t)} 
Js 

for some constants K and 5o. Moreover, (p'^{S) = mIS also satisfies 

R rSo 
a — 1 Js 



s 



(6.71) 



(6.72) 



as one easily sees by integration, and consistent with the fact that (j) solves (|6.7U|) 
with e = 0. Thus 



a 



- 1 Js 

13 . . rSo 



a — 1 



So""+^S" + e,5" / t-^dt + KS" 



/3 



a — 1 



1 — a 



Q—a-\-l na q 



(6.73) 



where we have applied the mean value theorem for integrals with e* = £(5"*) for 
some 5,, G (5, Sq). Now if we choose Sq = S"" for some < a < 1, then every term 
in ()6.73|) is higher order than S. That is, 



iS)< 



/3 + e, 



a — 1 

from which we conclude that 

<^2(5)-02(5) 



a — I 



S, 



lim 



0, 



as claimed. This completes the proof of Lemma IHl and thus Lemma [S] as well. □ 



6.4 The asymptotics for S > 1 

Equation (|6.18() implies that v = ^ = = patS = l. Thus, the TOV metric 
inside the Black Hole continues to the empty space Schwarzschild metric at n = 0, 
f = 0, S" = 1, an event horizon for the outer TOV metric in light of the fact 
that at 5 = 1, = ^ = 1, c.f. (I^TTj) . It follows that if the FRW density p 
is small at 5 = 1, then it makes sense in cosmology to approximate the FRW 
solution with p = for all times S > 1. Assuming this. Theorem El implies that 
for each < a < 1, the unique solution Uf^{S) continues to the zero pressure, 
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A; = 0, OS solution at 5" = = 1, the moment when the shock wave hes exactly 
one Hubble length from the FRW center. That is, at 5" = 1, the shock wave 
emerges from the White Hole event horizon of an ambient Schwarzschild metric 
as an outward propagating contact discontinuity that bounds a finite FRW mass, 
this being exactly equal to the total mass of the ambiant Schwarzschild metric 
into which it propagates. Thereafter the interface continues out to infinity along 
a geodesic of the Schwarzschild metric outside the Black Hole. 

We conclude that the OS solution gives the large time asymptotics of this 
new class of shock wave solutions that evolve inside of a Black Hole — and thus 
the explosion that begins at the Big Bang eventually settles down to a localized 
expansion that looks something like a giant supernovae, but on an enormous scale. 



7 Estimates for the Shock Position 

In this section we use the invariant region Q", (c.f. (|6.32|) ). to estimate the shock 
position fo at present time in terms of its position at the instant of the Big Bang 
S = 1/N = 0. Using this, we finish by estimating the time at which the shock 
wave first leaves the Black Hole in terms of the time at which the shock wave first 
becomes visible to an observer at the FRW center. 

Since the physical shock position at 5 = is f = 0, we begin by estimating 
fo = ro, R{to) = 1, in terms of r*, the value of the radial FRW coordinate at 
the instant of the Big Bang. In particular, the analysis shows that r^{S) has a 
limit = r^{a) as S — > 0, for any solution {r„{S),Ua{S)) of system (|5.,S5|) . (|5.36j) . 
So assume the FRW solution for the equation of state p = ap, < a < 1/3, is 
given, together with the shock trajectory Ua{S), that solves equation (|5.35j) . There 
then remains the one equation 1)5. 36() . leaving one free initial condition to impose. 
The next lemma gives an equivalent form of equation (|5.36|) in terms of the FRW 
variable r. 

Lemma 7 Equation iS.clb]) is equivalent to 

dr a — u r 



dS (l + fT)(l + 3u)5' 

for allO< S <1. 

Proof: Starting with (|5.3()j) . we can write 

df df dN 1 f 



(7.1) 



dS dN dS (l + 3u)S' 



(7.2) 



54 



Also, 



dr d(^) 1 df r dR , 
- - (7.3) 



dS dS RdS R'^dS' 
so using (|7.2|) we have 



dr / 1 \ r r dR 
dS " \l + 3u) S ~ R^' 

So now we need only compute dR/dS. For this, we use ()2.42|) . ()2.43|) to obtain 

3(1+ct) 

H = HqR —, 

and from (|5.1(ij) . 

ifV^ = N =1/S. 

These imply 



R=[SHy)''''-^' . (7.4) 
Using (|7.4|) together with (|7.2|) gives 



dR _ (1 + u) 
Is ~ (l + f7)(l + 3n) 



2 -2 

r 



—l (2+3a) 
3{1+<t) ^-3(1+^) 



(7.5) 



From (|7.4|) we get 



2=2 



i?3(l + (7) 



^0' ^ 

and using this in ()7.5() gives 



di? _ (1 + n) R 
dS (l + a)(l + 3n)5' 

Use this in (f7^ to get (f7l|) . □ 

We now use ()7.1|) to estimate tq in terms of r,,, assuming that the shock position 
^0 = '^o = V^/Hq at present time lies beyond one Hubble length. Then A^q > 1, 
and we have 
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1 



l/No = < 1, (7.7) 



c.f. (|2.7|) . (|5.17|) . Then integrating equation (|7.1j) we obtain, 

ro = r^e-'o ) {i+a)s''° ^ (7.8) 

where u denotes the function of S given by the trajectory u = Ua{S)) of (|5.H5j) . 
no-(O) = fj < 1/3, U(j(l) = 0, u'^{S) < 0, < 5 < 1. We can now use the invariant 
region to estimate a — Ua{S) for < 5 < 1. That is, the condition that the 
orbit Ua{S) lies in imphes that 

QaMS)) = a' 7^^^ <S< Min{l, hMS)} , (7.9) 
(1 + u)^ 

holds all along the orbit, where 

= {<r-n){l-3u) ^ {a-u){l-3u) ^ a-u 
^ ' 6u(l + n) - 6m - 6u ' ^ ' 

and we can take, c.f. (|6.33|) . 

„^ = ±. ,7.n, 

Now using (|7.11j) in the lower bound in (|7.9j) leads to 



1 + 3u - V 1 + 3n 
and applying this in (|7.8|) gives 



< ("I+Ji ) V3^V5 < V3^V5, (7.12) 



ro < r,e~^^°. (7.13) 



We now apply the upper bound in (|7.9jl . Note first that 



Min{l,h(uJS)} < Minil,^^-^^^^^-^} < Minil,^-^} , (7.14) 

L 6u J [ 6u ) 

where u is the (smallest) value of u at which h(u) = 1. (We need u, the smallest 
value of u that puts /^(m) < 1, to bound the factor u in the denominator of h{u).) 
We can estimate u as follows. First, note that by (|7.1()|) . h{u) = 1 is equivalent to 
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u = - — — — < 



7 + 3cr 7 + 3o- - 7 + So- 
so using the latter inequality to estimate the second term after the equality, we 
obtain 



u > 



a '6 f a •^'^ 



7 + 3a 7 + 3a \ 7 + 3a 



Using H7.15() in the second inequality in ()7.14|) . (i.e., ignoring for the moment the 
factor (1 — 3u) in the middle term of (|7.14j) ). implies that 

5<Mm(l,^^l <Mm|l,-^^l, (7.16) 
[ 6u ) [2a) 

which implies that 



a-u>^aS (7.17) 
holds all along the orbit u = UaiS). Using (|7.17j) in (|7.8|) gives the inequality 



fo>r,e3'"^«. (7.18) 

In the case a = 1 /3, we can improve the estimate ()7.18() by using the first inequality 
in (f7?Tl|l to obtain 

S < Afzn |l, < Afzn |l, , (7.19) 



18^ - 2 a 



and 



a-u> V2^, (7.20) 

in place of (|7.16j) and (|7.17|1 . Using (|7.20j) in (|7.8j) gives the improved inequality 
valid for cr = 1/3, 

ro > ne^^. (7.21) 

Putting H7.13() and H7.18() together, we obtain the following bounds for the shock 
position ro = fo in terms of the initial position r* that apply for < cr < 1/3 : 
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r,e^'"^« < ro < r^e^^^. (7.22) 
In the case p = ^p, we obtain the improved bounds 

r^e^^ < ro < r*et^. (7.23) 

The following Corollary follows immediately from H7.22() and ()7.23p in the case 
So = l: 

Corollary 2 Let f = fcrit denote the FRW shock position at the instant S = 1 
when the shock wave emerges from the Black Hole. Then for < a < 1/3 we have 

ne^" < fcrit < ne^, (7.24) 



while if a = 1/3 we have, 



^ - 2 
ne 4 < rcrit < ne2, (7.25) 



where r^, is the FRW radial coordinate of the shock wave at the instant of the Big 
Bang. 



Note that because 



S=^ = 5^. (7.26) 



multiplying 1)7.22(1 through by Hq gives 



—E, <^/S^<-E, <-E—, (7.27) 

mr^ Hor^ H^r^ 

and so we could use (|7.27() in (|7.22|) . (|7.23|) to obtain estimates involving Hq in 
place of Sq. In particular we have 



n < ro < r=,e(i+-'«o'-* , (7.28) 

which reproduces the OS result ro = r^ in the limit o" ^ 0. 

The final theorem gives an estimate for the number of Hubble lengths to the 
shock wave at the instant when it first becomes visible at the FRW center, as well 
as an estimate for the time it takes the shock wave to emerge from the Black Hole 
after it first becomes visible at the FRW center. 
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Theorem 9 Let r* = lim^^^o '"o-('S') denote the FRW position of the shock wave at 
the instant of the Big Bang, and assume < a < 1/3. Then the shock wave will 
first become visible at the center f = of the FRW spacetime at FRW time t = to, 
at the moment when the Hubble constant Hq = //(to) satisfies 



Hor, = ——, (7.29) 

i + 0(7 

( assuming R = 1 at t = to), and, the number of Hubble lengths \/Nq from the 
FRW center to the shock wave at time t = to satisfies 



2 ./3^(ijtM) 



1< < V^^o < e''^''^"^. (7.30) 

Furthermore, the time tcrit > to o^t which the shock wave will emerge from the 
Black Hole given that it first becomes visible at t = to is estimated by 



2 1^ trrii 2 2v^3f7 

e^'^ < — < — e i+- , (7.31) 



1 + 3cr ~ to ~ 1 + 3cr 
and by the better estimate 

e 4 < < 62, (7.32) 

^0 

in the case a = 1/3. 

Note that, for example, H7.3U() . H7.31() imply that at the OS hmit a = 0, 

0Vo = 2, hiii = 2, 
to 

and in the hmit a = 1/3, 

1.8<^<4.5, l<y^<A.5. 
to 

Note also that (|7.3fl|) and (|7.3H) imply that the shock wave will still lie beyond one 
Hubble length at the time t = to when it first becomes visible at the FRW center. 

Proof of Theorem Equation (|2.45|) implies that if the shock wave is first 
visible at t = to, -Ro = 1) then 
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which imphes (|7.29() . To verify H7.3U() . multiply equation H7.28() through by Hq, 
and use (f7^ and (|5T7|) . 

To verify (|7.31l) and (|7.32l) . let fcrit = TcntRcrit denote the shock position at 
S = N = \^ (the instant when the shock wave emerges from the Black Hole), and 
use and (f7?^ to write 



^ = = Hof^ra. (7.34) 

Then multiplying (|7.24|) and (|7.25|) through by Hq and using (|7.33|) in (|7.34|) gives 
()7.31|) and (|7.32|) . respectively. This completes the proof of Theorem El ^ 



8 Concluding Remarks 

We have constructed global exact solutions of the Einstein equations in which 
the expanding FRW universe extends out to a shock wave that lies arbitrarily 
far beyond the Hubble length. The critical OS solution inside the Black Hole is 
obtained in the limit of zero pressure, but the shock wave solutions have qualitative 
differences from the OS solution. For example, the shock surface f{t) tends to zero 
as t ^ 0, and the mass function at the shock M(f(t)) is finite for all t > 0, 
but unlike the OS solution, in the shock wave solution, M(f(t)) tends to infinity 
as t ^ 0. That is, the mass function is infinite at the instant of the Big Bang, 
but immediately becomes a finite decreasing function of FRW time, for all future 
times t > 0. Moreover, when p ^ 0, the directional orientation of the shock wave 
motion relative to the various observers is determined by the entropy condition- 
the entropy condition chooses the explosion over the implosion. For the entropy 
condition we take the condition that the pressure and density be larger behind the 
shock wave; that is, larger on the side that receives the mass flux. This condition 
implies that the shock is compressive, and is sufficient to rule out expansion shocks 
in classical gas dynamics, [T2] . 

One can ask the question, what is the solution like beyond the shock wave 
at any fixed instant of time inside the Black Hole? The answer is that the TOV 
energy density p{f) and the TOV mass M(f) as well, are both constant at each fixed 
"time" in the TOV metric beyond the shock wave, because f is the timelike TOV 
coordinate inside the Black Hole. This is no contradiction because what we identify 
(via shock matching) as the total mass function on the TOV side of the shock 
comes from the df^ component of the metric, which is timelike inside the Black 
Hole. However, this TOV "total mass" matches the FRW total mass continuously 
at the shock surface, and the FRW total mass has the physical interpretation 
as M = ^pf^, the integral of the energy density at each fixed time t in the 
FRW coordinates {t, r). Thus the evolution of the "total mass "is interesting and 
surprising inside the Black Hole. 

^^To see how M can be constant on the TOV side when it measures a total mass on 
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Throughout its expansion, the strength of the shock is on the order of the 
energy density p on the FRW side of the shock. At the moment when the shock 
wave Hes at the critical distance of exactly one Hubble length from the FRW 
center, the TOV density and pressure are zero, and thus we argue that if the FRW 
density is small as well, then the shock wave continues, (with small errors), to a 
zero pressure OS interface leaving the Black Hole at that instant. Thus the OS 
solution provides the large time asymptotics of these shock wave models. That is, 
the interface that marks the boundary of the FRW expansion continues out through 
the White Hole event horizon of an ambient Schwarzschild metric at the instant 
when the shock wave is exactly one Hubble length from the FRW center f = 0, 
and it then continues on out to infinity along a geodesic of this Schwarzschild 
metric outside the Black Hole. These solutions thus indicate a scenario for the 
Big Bang in which the expanding universe emerges from an explosion emanating 
from the White Hole singularity inside the event horizon of an asumptotically fiat 
Schwarzschild spacetime of finite mass. The model does not require the physically 
implausible assumption that the uniformly expanding portion of the universe is 
of infinite mass and extent at every fixed time, and it has the nice feature that 
it embeds the Big Bang singularity of cosmology within a larger spacetime, the 
Schwarzschild spacetime. Moreover, the model also allows for arbitrarily large 
densities to exist over arbitrary numbers of Hubble lengths early on in the Big 
Bang, a prerequisite for the standard physics of the Big Bang at early times. 

One might ask how an observer near the FRW center would first detect evidence 
of such a cosmic shock wave. Since the shock wave emerges from the Big Bang 
beyond the Hubble length, the model would imply a uniform expansion throughout 
a region that is initially well beyond the light cone of an observer positioned 
near the FRW center. If the shock wave were initially far enough out, then the 
uncoupling of matter from radiation at about 300, 000 years after the Big Bang 
would produce an extended region with a uniform background radiation field. 
This would persist until roughly the time when the Hubble length catches up to 
the shock wave, a time determined by the initial conditions. The influence of the 
solution beyond the shock wave would propagate into this radiation field at the 
speed of light, first appearing to an observer that is off center on the FRW side 
of the shock as a disturbance in the background radiation field at a point in the 
sky in the direction nearest the shock wave, and the disturbance would grow from 
that time onward. 

These exact shock wave solutions give the global dynamics of strong gravita- 

thc FRW side, consider shock matching in {t.f) coordinates. Then since the FRW mass 
M depends only on t, while the TOV mass depends only on f, one of them depends on 
the timelike coordinate and one on the spacelike coordinate in the {t, f) coordinates at the 
shock. Thus, in the Einstein equations, the nonzero M derivative ends up in a different 
equation on each side of the shock, giving the physical total mass on the FRW side. But 
on the TOV side, the derivative M' is equated to the pressure, and thus doesn't have the 
same interpretation as an integral of the energy density. 
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tional fields in an exact solution, the dynamics is qualitatively different from the 
dynamics of solutions when the pressure p = 0, and the solution suggests a Big 
Bang cosmological model in which the expanding universe is bounded throughout 
its expansion. Surprisingly, unlike shock matching outside the Black Hole, the 
equation of state p = of early Big Bang physics, plays a special role in the 
equations, and for this equation of state alone, the behavior of the shock wave at 
the instant of the Big Bang is distinguished. But these solutions are only rough 
qualitative models because the equation of state on the TOV side is determined 
by the equations, and therefore cannot be imposed. That is, the TOV density p 
and pressure p only satisfy the loose physical bounds < p < p ; and on the FRW 
side, the equation of state is taken to be p = ap, a = const., < a < 1. We take 
these bounds as implying that the equations of state are qualitatively reasonable. 
The entropy condition, p > p, p > p, (that the density and pressure be larger on 
the side that receives the mass flux), implies that the shock wave is compressive, 
and this fixes a time orientation for solutions, and determines a unique solution. 

However, we expect that these solutions will capture the gross dynamics arising 
when more general equations of state are imposed. In fact, we suggest that the 
global dynamics described in these solutions could only be discovered within a class 
of exact solutions in which simplifying assumptions are made. For more general 
equations of state, other waves, (e.g. rarefaction waves), would need to be present 
to meet the conservation constraint, and thereby mediate the transition across the 
shock wave. Such transitional waves would be pretty much impossible to model in 
an exact solution. 

Finally, we note that because Einstein's theory by itself does not choose an 
orientation for time, it follows that if we believe that a Black Hole can exist in the 
forward time collapse of a mass through an event horizon as t — > oo, (the time t as 
observed in the far field), then we must also allow for the possibility of the time 
reversal of this process, a White Hole explosion of matter through an event horizon 
coming from t — oo. These solutions might be relevant in explaining astropysical 
systems, such as galaxies and stellar associations whose expansions appear so 
great as to have emerged from an event horizon at earlier times-an impossibility 
if one only allows the time orientation of a collapsing Black Hole, and not its 
time reversal, the expanding White Hole. Of course, this naturally leads one to 
wonder if there is a connection between the mass that mysteriously disappears 
into Black Hole singularities, and the mass that mysteriously emerges from White 

"'^^The time orientation of a solution must be selected based on an extra condition, such 
as an entropy condition for shocks, because the Einstein equations and the compressible 
Euler equations, taken by themselves, are both time reversible, j7l I12j . Thus the entropy 
condition for the shock is what determines the time orientation for the global dynamics 
of the solutions we construct: the FRW metric expanding outward behind a shock wave 
emanating from a White Hole is entropy satisfying, while its time reversal, the FRW metric 
contracting into a Black hole, is entropy violating. We find it interesting that the entropy 
condition determines a unique solution in the large. 
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Hold singularities. 
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